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13.  ABSTRACT 


Three  computer  programs  are  described  which  allow  the  study  of  optical 
properties  of  two  dimensional  laser  resonators.  Two  of  the  programs 
treat  the  propagation  of  an  electric  field  inside  an  optical  resonator  (there 
arc  two  methods  for  propagating  an  arbitrary  electric  field),  and  couple 
the  propagating  fields  to  a  specified  laser  medium  which  may  be  nonuni¬ 
form,  transversely  flowing,  saturable,  and  inhomogeneous.  The  propaga¬ 
tion  technique  includes  diffraction  effects  throughout  the  laser  medium. 

Both  programs  provide  output  beam  characteristics  and  output  power.  The 
third  program  calculates  far  field  beam  characteristics  •  jr  the  laser  output. 


Some  of  the  features  offered  the  user  are  a  restart  facility  for  long  runs, 
arbitrary  location  of  the  medium  in  the  resonator  cavity,  and  capacity  for 
mirror  tilt  on  either  or  both  mirrors. 
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FOREWORD 


This  manual  describes  three  computer  programs  designed  to 
study  the  optical  properties  of  laser  resonators.  The  first  two  are 
both  Step-by-Step  Computer  Programs,  termed  SSCP  and  SSCP(ECS), 
where  the  latter  is  based  on  an  Expanding  Coordinate  System.  Both 
programs  simulate  the  propagation  of  one  or  more  waves  between 
mirrors  in  a  laser  cavity,  differing  primarily  in  the  coordinate  system 
used  to  describe  the  propagation.  Experience  has  shown  the  SSCP(ECS) 
program  to  be  more  efficient  in  most  cases,  requiring  less  computer 
time  and  memory  to  solve  equivalent  problems,  and  for  this  reason 
only  the  SSCP(ECS)  is  fully  documented.  The  SSCP  package  is  docu¬ 
mented  only  in  Section  I  and  Appendix  C. 

The  third  program  is  designed  to  take  an  output  waveform  from 
the  laser  cavity  (computed  by  one  of  the  two  programs  above)  and 
propagate  it  to  the  far  field  to  determine  the  resulting  intensity  distri¬ 
bution.  Appendix  B  contains  a  complete  description  of  the  program, 
including  its  mathematical  basis  anu  a  description  of  the  required 
input. 

This  report  represents  the  second  of  two  user's  manuals  which 
have  been  prepared  under  Contract  DAAH01  -73-C-0298.  The  first  was 
a  user's  guide  for  laser  resonator  modeling  using  the  Fresnel  integral 
technique.  This  work  was  monitored  by  William  Martin  of  the  Propulsion 
Directorate. 
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Sy  mbol 


A  A 


A  SR 


D1A1,  D1A2 


Defined  by  eq.  (14) 

t  for  cylindrical  resonators,  4  for  square 
resonators,  O  for  2-dimensional  calculations 

Speed  of  light 

Diameter  of  mirror  1  (output  mirror)  or 
mirror  2 


DXMIR  1 ,  DXMIR2 

G. 

J 

“l 

«2 

GZERO 


Mesh  spacing  at  mirror  1  or  2  for  x 
coordinate 

th 

Small-signal  gain  at  j  mesh  point 

Resonator  parameters  defined  by  eq,  (16) 

Resonator  parameter  defined  by  eq.  (16) 

Small-signal  gain  for  analytical  gain 
subroutine 


Local  flux  density  for  two-dimensional 
resonator 


Local  flux  density  for  pseudo-three- 
dimensional  resonator 


K  Wave  number  (2t r/\) 

o 

M  Geometric  magnification  of  unstable 

resonator 


NX 


Number  of  propagating  field  points  in 
x  direction 


NZ 


Number  of  propagation  steps  between 
mirrors 


OMEGA  1 ,  OMEGA2  1/e  truncation  distance  for  mirror  edges 

RCURV  1 ,  RCURV2  Radius  of  curvature  for  mirror  1  and 

mirror  2 
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Sy  mbol 


REFL1,  REFL2 
T  LOSS  1 

y 

A 

u 

U  m 
J 

U2d 

VR1,  VR2 
x 

XIZERO 

X  WIDTH 

ZLEN 

ZLENMD 

ZMDLOC 


Greek 


AP 


2d 


Ax 


Az 

An 

An 


6x 

X 


Mirror  reflectivity 
Transmission  loss  for  mirror  1 
Electric  field  vector 
Phasor  for  electric  field  U 

Free-space  electric  field  at  x  =  jAx  and  z  =  mAz 

Electric  field  in  two-dimensional 
resonator 

Defined  in  eq.  (15) 

Transverse  coordinate  for  electric  field 

Saturation  parameter  for  analytical  gain 
subroutine 

Extent  of  x  dimension  at  start  of  each 
roundtrip 

Distance  between  mirrors 
Gain  medium  length 

Distance  from  mirror  1  to  center  of  gain 
medium 

n2-  1 

Power  contained  in  strip  of  Ax  width 
Incremental  step  size  in  x  coordinate 
Incremental  step  size  in  z  coordinate 

n  -  1 

th 

nj  -  1 ,  where  n  =  refractive  index  at  j 
field  point 

Defined  in  eq.  (7) 

Wavelength 


vm 


Laplacian  operator  in  rectangular 
coordinates 

Laplacian  operator  for  x  coordinate  only 

Angular  coordinate  for  cylindrical 
resonators 


Optical  frequency 


I.  SYSTEM  DESCRIPTION 

A .  Introduction 

The  purpose  of  the  Step-by-Step  Computer  Program  (SSCP)  is 
to  study  the  optical  properties  of  laser  resonators.  The  optical  prop¬ 
erties  of  a  laser  system  are  dependent  upon  many  factors  which  are 
interrelated  in  a  complicated  manner.  This  makes  it  impossible  for  an 
experimenter  to  determine  the  best  choice  of  medium  and  resonator 
parameters  for  producing  an  optimum  output  beam  without  an  expen¬ 
sive  experimental  trial  and  error  analysis  or  without  conducting  a 
realistic  analytical  study.  The  SSCP  is  useful  to  the  experimenter 
because  it  can  greatly  narrow  the  range  of  parameters  from  which  he 
has  to  choose.  By  modeling  the  laser  medium  and  the  optical  resona¬ 
tor  for  the  device  under  analysis,  he  can  vary  the  model  parameters 
which  represent  physical  variables  and  choose  a  range  of  parameters 
that  give  the  best  beam  quality. 

The  experimenter  has  the  choice  of  making  his  analytical 
analysis  using  the  Step-by-Step  Computer  Program  (SSCP)  described  in 
this  user's  manual  or  the  Fresnel  Integral  Computer  Program  (FICP) 
described  in  the  use/'s  manual  for  the  Two-Dimensional  Fresnel 
Integral  Computer  Program.  The  choice  is  dependent  upon  the  size 
of  the  resonator's  Fresnel  Number  and  the  homogeneity  of  the  laser 
medium.  The  SSCP  must  be  used  for  laser  resonators  with  Fresnel 
numbers  larger  than  sixty,  for  a  transverse  flowing  medium,  when  the 
medium  is  highly  gain  saturated,  or  when  significant  density  variations 
exist. 

The  SSCP  is  a  very  versatile  tool  for  the  experimenter  to  use 
in  the  analysis  of  large  Fresnel  number,  inhomogeneous  laser  devices. 
Once  he  has  optimized  the  nonoptical  parameters  to  obtain  the  largest 


Two-Dimensional  Fresnel  Integral  Compute  :  Program,  User's 
Manual,  Hughes  Research  Laboratories,  Contract  NoT  DAAH01-73 
C-0398,  May  1973. 
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amount  of  energy  in  the  medium,  he  must  select  an  optical  resonator 
which  will  optimize  the  far  field  intensity  distribution  (in  most 
instances  this  means  obtaining  the  largest  far  field  flux  density).  The 
SSCP  will  allow  the  experimenter  to  couple  up  to  6  laser  waves  to  a 
particular  laser  medium  (except  for  empty  cavity  calculations  the 
SSCP  must  be  furnished  with  a  kinetics  subroutine).  The  medium  sub¬ 
routine  can  represent  a  single  line  laser  system  such  as  a  CO  laser  or 
a  multiline  laser  such  as  an  HF  laser  and  it  must  supply  saturated  gain 
and  refractive  index  for  each  local  propagating  field  point.  It  will 
accurately  account  for  diffraction  of  each  propagating  wave  as  they 
propagate  through  the  laser  medium.  By  examining  each  wave  and 
their  composite  (obtained  from  their  averaged  intensities)  the  user  can 
determine  if  the  laser  will  be  single  or  multimode,  the  total  roundtrip 
loss  or  gain  for  each  wave,  the  output  power  of  the  laser,  the  intensity 
and  phase  distribution  for  each  wave  in  both  the  near  and  far  field,  and 
an  averaged  intensity  distribution  for  both  the  near  and  far  field  if  the 
laser  is  operating  multimode.  All  the  above  information  „  ,\n  be  both 
printed  and  plotted  after  each  roundtrip  pass  through  the  resonator  (the 
operator  has  the  option  to  print  and  plot  after  any  given  roundtrip 
pass).  Furthermore,  the  user  has  the  option  of  including  mirror  tilt, 
selecting  the  axial  location  of  the  medium  inside  the  resonator,  and 
selecting  mirror  reflectivity  and  mirror  transmission  loss. 

The  accuracy  of  the  optimization  results  using  the  SSCP  depends 
upon  how  realistically  the  medium  and  resonator  models  describe  the 
actual  laser  device.  Since  any  model  must  make  certain  approxima¬ 
tions,  the  optical  properties  as  predicted  by  the  SSCP  will  never  be  in 
exact  agreement  with  measured  results.  Therefore,  it  is  important  for 
anyone  using  the  SSCP  to  understand  the  techniques  used  in  the  analysis 
and  its  limitations.  In  the  remainder  of  this  section  the  techniques  used 
to  analyze  laser  resonators,  the  type  of  information  available  from  the 
SSCP,  the  limitations  of  the  computer  program,  and  the  factors  govern¬ 
ing  input  data  to  the  program  will  be  described. 
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ihe  approach  that  is  taken  to  solve  for  the  transverse  modes  of 
a  laser  resonator  is  one  which  closely  resembles  the  actual  develop¬ 
ment  of  the  modes  from  the  laser's  spontaneous  emission.  An  arbi¬ 
trary  input  wave  (the  SSCP  has  a  choice  of  a  wave  with  uniform  or 
gaussian  amplitude  distribution)  is  propagated  back  and  forth  between 
two  resonator  mirrors.  As  the  wave  traverses  between  the  mirrors  it 
interacts  with  the  laser  medium.  As  it  interacts,  its  amplitude  and 
phase  are  adjusted  according  to  the  gain  and  refractive  index  of  the 
medium.  This  iterative  procedure  is  repeated  until  a  stable  solution  is 
obtained.  The  wave  represented  by  the  stable  solution  is  the  steady 
state  transverse  mode  or  modes  of  the  laser.  A  measure  of  the  beam 
quality  is  obtained  by  transforming  the  steady  state  solution  (near  field 
of  the  laser)  into  the  far  field  and  comparing  it  with  the  resulting  far- 
field  distribution  for  a  uniform  near-field  distribution. 

1 .  Approximate  Solution  to  the  Wave  Equation 

Confining  our  attention  to  the  case  in  which  orthogonal 
field  polarizations  are  not  mixed  by  the  laser  medium,  the  electric 
field  vector  U  within  the  resonator  may  be  represented  by 

U(x,  z,  t)  =  i  Re[U(x.  z)  exp  (iKQZ  -  iwt)](  (1) 

where  i  is  a  unit  polarization  vector,  Kq  is  the  free-space  wave  num¬ 
ber  (Kq  =  2 tt7M>  w  is  the  optical  frequency  in  rad/sec  (  w  =  cKq),  and 
the  phasor  U  is  a  complex  number  depending  on  the  position  coordi- 

A 

nates  x,  z.  In  a  laser  resonator,  U(x,  z)  satisfies  the  wave  equation. 

(VZ  +  KQ2n2)  U(x,  z)  =  0  (2) 

where  V2  is  the  Laplacian  in  rectangular  coordinates,  and  n  is  the 
refractive  index  (contains  both  real  and  imaginary  terms).  By  con¬ 
verting  the  above  differential  equation  into  an  integral  equation  and 


3 


using  the  Green's  function  for  the  Helmholtz  equation,  we  can  obtain  the 
following  equation: 

U(x,  z  +  Az)  =  Uq(x,  z  +  Az)  exp  ^i  K^AnAz  +  0(cAz^/p^)j  (3) 

A 

where  U  (x,  z  +  Az)  is  the  resulting  electric  field  after  propagating 

a  °  i 

U(x,  z)  through  free-space  for  a  distance  Az,  t  =  n  -  1,  p  is  the 

shortest  characteristic  length  in  the  z  direction  over  which  the  field 

Z  Z 

amplitude  or  phase  changes,  and  0(tAz  /p  )  is  a  truncation  error 

A 

term.  In  other  words,  if  we  know  both  the  electric  field  U  and  the  local 

refractive  index  of  the  medium  (assumed  to  be  constant  for  a  distance  Az) 

in  the  z  plane,  and  if  we  know  what  the  free-space  electric  field  U  is  in 

o 

A 

the  z  +  Az  plane,  we  can  determine  the  electric  field  U  in  the  z  +  Az 
plane. 

In  a  laser  resonator  model,  the  medium  as  well  as  the  initial 
propagating  field  is  assumed  to  be  known.  Therefore,  what  remains  to 
be  solved  is  how  to  determine  the  free-space  electric  field  at  z  +  Az 

A 

from  a  known  electric  field  U  at  z. 

Z.  Propagation  Algorithm  (Two-Dimensional) 

The  two-dimensional  free-space  wave  equation  assum¬ 
ing  constant  Kq  is  written  as 


(i2Ko-|+  Vt2)6o<x'2>  =  0  <4> 

z  z  z 

where  =  3  / 3x  .  Using  the  following  difference  equations  of  DuFort 
and  Frankel  (Ref.  Z) 


R.  D.  Richtmyer  and  K.  W.  Morton,  Difference  Methods  for  Initial- 
Value  Problems,  Interscience  Publishers,  Second  Edition,  1967. 
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u 

<§W 


* 


I 


A 

au 


_ o  _  1  m+1  ..  m-  1  \ 

dz  ~  2Az  \Uj  "  Uj  / 


(5a) 


2  A 

a  u 


ax  Ax 


-Mu."'.  -  u  m+1  -  u  m-‘  +  u  m,) 

Av^  '  J+1  J  J  J-1/ 


(5b) 


where 


x  =  jAx 
z  =  mAz 


1  <  j  <  NX 
1  <  m  <  NZ 


(6) 


[Eq.  (5)  includes  both  the  real  and  imaginary  terms  of  U  .]  We 
obtain  after  substituting  eq.  (5)  into  eq.  (4)  the  following  propagation 
algorithm: 


U  m+1  =  — M  (b  m-‘  +  itx(un;,  -  2U  n'-‘  +  5”  ) 
J  1  +  tx2  L 1  \  J+1  >  >-d 

,  2  / *  m  A  m-1  ,  *m  \  \ 

+  „x(Uj+l  -U.  +U..JJ 


+  Truncation  Terms 


(7) 


where  ftx  =  Az/ ^Kq(Ax)^ j .  Equation  (7)  allows  the  field  at  z  +  Az  to 

be  computed  from  known  field  values  at  z  and  z  -  Az. 

The  difference  method  of  DuFort  and  Frankel  is  explicit  and  is 

unconditionally  stable  (solution  always  converges  independent  of  ox, 

and  Az)  in  unbounded  regimes.  From  experience  and  an  examination  of 

2 

the  truncation  terms,  it  has  been  determined  that  if  Az  £  0.  5Kq(Ax)  , 
and  Ax  and  Az  are  less  than  the  smallest  characteristic  length  in  the  x 
and  z  directions,  respectively,  eq.  (7)  will  accurately  propagate  an 
arbitrary  field  through  free  space  with  a  gain  or  loss  of  less  than  1  or 
2%  for  each  roundtrip  pass  (the  user  should  note  that  these  errors  do 
not  accumulate  from  pass  to  pass). 
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3. 


Calculation  Procedure 


f  i 


The  calculation  technique  for  analyzing  laser  resonators 
using  the  step-by-step  propagation  algorithm  given  in  eq.  (7)  is  very 
similar  to  the  technique  described  in  Ref.  1;  the  main  difference  is  the 
propagation  algorithm  [Ref.  1  uses  a  Fresnel  integral  to  propagate  the 
electric  field  back  and  forth  in  the  resonator,  while  here  propagation 
is  accomplished  using  the  numerical  difference  technique  in  eq.  (7)]. 

The  approach  taken  is  to  consider  a  propagating  wave  or  waves 
reflecting  back  and  forth  between  two  cavity  mirrors.  It  will  be 
assumed  for  the  moment  that  there  is  no  laser  medium  present.  (This 
condition  gives  the  empty  cavity  or  nonperturbed  mode.)  The  follow¬ 
ing  calculation  sequence  is  used: 

Empty  Cavity 

1.  Reflect  input  wave  or  waves  (a  maximum  of  six 
waves  can  be  used)  from  mirror  1.  This  step 
corrects  the  field  amplitude  for  mirror  reflectiv¬ 
ity  and  mirror  size  and  corrects  the  phase  of  the 
field  for  mirror  curvature  and  angular  mirror 
tilt.  (For  perfectly  aligned  resonators,  mirror 
tilt  is  zero. ) 

2.  Propagate  one  step  to  z  =  Az  using  the  following 
initial  propagation  algorithm 

A  1  *0  /a  o  a  q  a  o\ 

U.  =U  +i6x  U.  ,  -  2U.°  +  U...  )  (8) 

J  J  V  J"1  )  J+1  / 

3.  Propagate  to  mirror  2  (see  Fig.  1)  using  eq.  (7) 
repeatedly  for  NZ  -  1  steps. 

4.  Reflect  from  mirror  2.  This  step  corrects  the 
amplitude  and  phase  of  the  field  for  the  effects  of 
mirror  2. 

5.  Propagate  one  step  using  eq.  (8). 

6.  Propagate  to  mirror  1  using  eq.  (7)  repeatedly 
for  (NZ  -  1)  steps. 
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7. 


A 

Store  real  and  imaginary  terms  of  U:  on  a 
RSTART  tape.  These  values  will  be  used  later 
in  a  FAR  FIELD  PROGRAM  to  compute  the  far 
field  intensity  distribution. 

8.  Renormalize  field  amplitude  at  mirror  1.  This 
step  multiplies  the  field  at  mirror  1  by  a  constant 
so  that  at  this  location  the  area  under  the  distri¬ 
bution  curve  remains  fixed  from  pass  to  pass. 

9.  Repeat  steps  1  through  8  for  NPASS  passes. 

By  following  steps  1  through  9,  the  many  modes  that  went  into  defining 
the  input  wave  will  reduce  to  a  single  mode,  and  this  mode  will  be  the 
lowest  loss  transverse  mode  of  the  resonator.  This  mode  is  also  the 
lowest  loss  mode  for  three-dimensional  square  mirror  geometry  and 
should  be  interpreted  as  applying  to  a  cross  section  through  a  resona¬ 
tor  having  square  mirrors  with  spherical  surfaces.  It  should  be  noted 
that  when  a  saturable  gain  medium  exists,  this  equivalence  is  no 
longer  true. 

Figure  2  shows  the  active  laser  medium  inserted  into  the  calcu¬ 
lation  procedure.  Unlike  the  FICP,  the  SSCP  is  not  restricted  to  low 
Fresnel  number  calculations  and  therefore  can  treat  diffraction  in  the 
gain  medium.  The  total  gain  length  of  the  medium  is  ZLENMD,  where 
ZLENMD  can  be  equal  to  or  less  than  ZLEN.  The  location  of  the  cen¬ 
ter  of  the  laser  medium  is  located  ZMDLOC  distance  from  the  mir¬ 
ror  1 . 

When  propagating  through  the  laser  medium,  the  phase  and 
amplitude  of  the  field  is  adjusted  after  each  propagation  step  according 
to  eq.  (3).  Rewriting  eq.  (3)  in  terms  of  the  local  power  gain  Gj  and 
the  local  index  of  refraction  An.,  results  in  the  following: 


U(jAx,  z  +  Az)  =  UQ(jAx,  z  +  Az)  exp 


[(-/ +  ‘  K<,Anj) 


(9) 


8 


2785-2 


Y 


Fig.  2.  Resonator  geometry  with  laser  medium. 


(G  can  be  positive  or  negative  depending  on  whether  the  medium  has 
gain  or  loss),  Az  is  the  step  size  in  the  gain  medium  and  is  equal  to 
ZLEN/NZ  and  Am  is  the  local  refractive  index  difference.  The  local 
power  gain  G;  can  be  supplied  by  a  separate  gain  subroutine  such  as 
the  Lockheed  2-D,  Di  -CO^  (Ref.  3)  mixing  program  or  by  a  built-in 
gain  subroutine  GAINA  which  assumes  the  gain  to  be  uniform  and  homo¬ 
geneously  broadened.  The  latter  gain  routine  is  expressed  as 


G. 

J 


GZERO 


1  + 


Ux(jAx)| 

XIZERO 


(10) 


where  GZERO  is  the  small-signal  gain,  XIZERO  is  the  saturation  flux 
density,  and  UT(jAx)  is  total  field  amplitude  (forward  plus  backward 
traveling  wave).  If  the  medium  is  assumed  to  have  uniform  density,  then 
Am  =  0.  Since  this  type  of  gain  medium  cannot  couple  laser  transitions, 

GZERO  and  XIZERO  are  the  same  for  each  propagating  wave. 

The  calculation  procedure  when  the  two-dimensional  laser 

medium  is  included  is  as  follows: 


Medium  Included: 


1.  Reflect  input  wave  or  waves  from  mirror  1. 

2.  Propagate  to  edge  of  laser  medium  using  eqs.  (7) 
and  (8). 

3.  Sum  present  field  intensity  with  field  intensity 
from  previous  pass  traveling  in  opposite  direc¬ 
tion  (stimulated  emission  depends  on  total  local 
intensity). 

4.  Call  gain  subroutine.  This  routine  furnishes 
gain  and  refractive  index  values  for  each  field 
point. 

5.  Propagate  one  step  in  the  z  direction  using 
eq.  (7). 

6.  Adjust  the  Uj  field  points  for  gain  and  density 
variations  according  to  eq.  (9). 


3 

Chemical  Laser  Analysis  Development  Program,  Two-Dimensional 
Fresnel  Integral  Laser  and  Mixing  Program  User's  Guide,  Vol.  II, 
Lockheed  Missile  and  Space  Company,  Contract  DAAH01 -73-C-0232, 
Technical  Report  RK-CR-73-2,  October  1973. 
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7.  Store  the  intensity  (|Uj|^)  for  each  field  point. 

(The  total  number  of  times  the  intensities  are 
stored  in  the  z  direction  is  limited  to  21.  If  the 
number  of  propagation  steps  in  the  medium 
exceeds  21,  the  program  automatically  skips 
storing  every  step  in  order  not  to  exceed  this 
number.  ) 

8.  Repeat  steps  5  and  6  for  total  of  IGSKIP  times. 

The  parameter  IGSKIP  requires  that  IGSKIP 
propagation  steps  be  made  before  the  gain  sub¬ 
routine  is  called. 

9.  Repeat  steps  3  through  8  until  far  side  of  laser 
medium  is  reached. 

10.  Propagate  to  mirror  2  using  eq.  (7). 

11.  Reflect  from  mirror  2. 

12.  Repeat  steps  2  through  9. 

13.  Propagate  to  mirror  1  using  eq.  (7). 

14.  Store  amplitude  and  phase  of  field  points  on  the 
RSTART  tape.  (This  data  is  used  in  the  FAR 
FIELD  PROGRAM  to  compute  the  far  field 
intensity  distributions.) 

15.  Repeat  steps  1  through  14  for  NPASS  passes. 

Laser  power  can  be  coupled  out  of  the  laser  at  either  or  both 
ends  of  the  resonator  if  the  beam  spills  past  the  edges  of  the  resonator 
mirrors.  The  power  lost  in  this  way  can  be  determined  by  taking  the 
beam  power  before  and  just  after  the  beam  hits  a  mirror  and  taking 
their  difference  (mirror  reflectivity  is  assumed  to  be  1).  Useful  laser 
power  is  assumed  to  occur  at  mirror  1  and  this  can  be  from  both 
spillage  and  transmission  through  the  mirror  (REFL1  <  1.).  It  is  this 
combined  beam  that  is  used  to  determine  the  far  field. 

4.  Mirror  Edge  Tapering 

The  truncation  error  terms  in  eq.  (7)  contain  the  fol¬ 
lowing  partial  derivatives. 


1 1 


Truncation  Error 


i 


Because  of  the  partial  derivatives  with  respect  to  x,  we  cannot  accu¬ 
rately  treat  infinitely  steep  field  gradients  that  occur  at  sharp  mirror 
edges;  therefore,  it  is  necessary  to  use  edge-tapering.  In  our  step- 
by-step  resonator  model,  the  mirror  reflectivity  at  the  edge  of  the 
mirror  is  given  by  the  following  gaussian  function: 


Mirror  Reflectivity  =  REFL  exp  ^-(x  -  XQ)^/OMEGA^j  ,  for  jx|  >  xq 

(1Z) 

where  REFL  is  the  power  reflectivity  of  the  mirror,  Zxq  is  the  mirror- 
width  describing  the  region  of  uniform  reflectivity,  x  =  o  at  mirror  cen¬ 
ter,  and  OMEGA  describes  the  distance  of  truncation.  (When  the  operator 
sets  DIA1  or  DIAZ,  the  program  automatically  computes  xq  such  that 

x  =  DIA/Z  -  v/o.  7  •  OMEGA.  This  fixes  the  effective  mirror  diameter 
o 

between  the  half  power  points  in  the  taper  region.  ) 

C.  Step-by-Step  Technique  (Expanding  Coordinates) 

The  discussion  up  to  this  point  has  outlined  the  technique  for 
solving  for  the  transverse  modes  of  a  laser  resonator  using  the  step- 
by-step  propagation  technique.  As  stated  in  the  introduction  the  step- 
by-step  technique  is  best  suited  for  large  Fresnel  number  waves. 

Thus,  the  step-by-step  model  should  be  best  suited  for  studying  large 
Fresnel  number  resonators.  However,  this  is  not  the  case.  Larger 
Fresnel  numbers  require  increasingly  larger  numbers  of  mesh  points 
and  therefore  more  computer  time.  This  is  due  to  the  beam  curvature 


1Z 


which  is  introduced  into  the  wave  every  time  the  beam  bounces  off  a 
resonator  mirror.  For  example,  in  a  confocal  unstable  resonator 
using  a  convex  and  concave  mirror  (laser  output  is  obtained  at  the 
convex  mirror),  the  beam  that  is  reflected  from  the  output  mirror  has 
a  radius  of  curvature  equal  to  that  of  the  mirror  and  thus  its  effective 
Fresnel  number  is  much  lower  than  the  Fresnel  number  of  the  resona¬ 
tor.  However,  by  using  an  expanding  coordinate  system,  we  can 
remove  the  dependence  of  the  mirror  curvature  on  the  number  of  mesh 
points  required  to  accurately  propagate  the  beam.  We  can  propagate 
the  beam  back  and  forth  inside  the  resonator,  assuming  it  is  nearly  a 
plane  wave,  and  then  by  using  a  suitable  coordinate  transformation 
determine  the  actual  beam  distribution  including  phase  information. 

The  SSCP  can  treat  diffraction  within  the  laser  medium,  but  it 
is  limited  to  studying  resonators  with  Fresnel  numbers  of  approxi¬ 
mately  60  or  less.  However,  by  going  to  an  expanding  coordinate  sys¬ 
tem  the  Step-by-Step  Computer  Program  (Expanding  Coordinate 
System)  SSCP(ECS)  does  not  have  an  upper  limit  on  the  Fresnel  num¬ 
ber.  Because  of  the  importance  of  large  Fresnel  numbers  in  high 
power  laser  systems,  this  user's  manual  will  be  mainly  devoted  to  the 
SSCP(ECS).  Features  which  are  explicitly  different  in  the  SSCP  will 
be  described  in  Appendix  C. 

1.  Propagation  Algorithm  (Two-Dimensional) 

The  propagation  algorithm  for  the  expanding  coordinate 
system  is  similar  to  the  algorithm  given  in  eq.  (7)  and  is  expressed  as 
follows. 


Chemical  Laser  Mode  Control  Program,  Semiannual  Technical 
Report,  Hughes  Research  Laboratories,  Contract  No.  DAAH01-73- 
C-0290,  13  November  1972  through  12  May  1973. 
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+  6x 

+  Truncation  Terms 


(13) 


where  AA  is  a  variable  which  is  a  function  of  z.  AA  is  defined  as 


AA  ”  ( 1  "  2(z  +  VR)  ) 


(14) 


where  VR  is  shown  in  Fig.  3  as  the  virtual  image  source  distance  with 
z  measured  from  the  mirror  as  the  beam  propagates  away  from  it. 


The  virtual  source  distance  V  l  is  a  function  of  the  physical  parameters 
of  the  resonator  and  is  given  as^: 


VR1.2  = 


1  -  (gjg2)  -  1  +(gi  ?) 


-1 


1,2' 


2  -(g/1 


ZLEN 


(15) 


w 


here  g,  are  the  resonator  g  parameters  defined  as 
1 9  £ 


,  ZLEN 

gl,  2  RCURV 1 ,  2  ' 


(16) 


Because  of  the  expanding  coordinate  system  6x  is  no  longer  a  constant 
since  Ax  is  a  function  of  z.  Ax  is  expressed  as  follows. 


A.  E.  Siegman,  "Unstable  Optical  Resonators  for  Laser  Applica¬ 
tions,  "  Proceedings  IEEE,  53,  277-287,  March  1965. 
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Ax  =  DXMIR1  |l  +  Y|^y)  propagating  from  mirror  1  to  mirror  2 

Ax  =  DXMIR2  ^1  +  YR2.~)  ProPa8at*nS  ^rom  mirror  2  to  mirror  1 

(17) 


where  DXMIR1  and  DXMIR2  are  the  values  for  Ax  at  mirror  1  and 
nirror  2,  respectively,  and  z  is  again  measured  from  the  mirror  that 
the  beam  is  propagating  away  from. 

As  a  result  of  the  coordinate  transformation  the  SSCP(ECS) 
differs  from  the  SSCP  in  that  beam  expansion  or  contraction  is  con¬ 
trolled  by  the  increase  or  decrease  of  Ax  in  eq.  (17)  and  not  by  the 
phase  curvature  placed  on  the  beam  upon  reflection  from  a  mirror  as 
is  the  case  for  the  SSCP.  In  the  SSCP(ECS),  the  transformed  mirrors 
are  flat  (infinite  radius  of  curvature).  The  operator,  when  inputting 
the  mirror  radii  of  curvature,  must  supply  the  physical  curvatures  to 
the  program  (these  are  used  to  compute  the  virtual  image  distances  in 
eq.  15),  but  the  program  treats  the  mirrors  as  if  they  were  flat  when 
the  beam  reflects  off  their  surfaces.  The  divergence  or  convergence 
that  would  be  caused  by  the  actual  curvature  of  the  mirrors  is  treated 
by  the  coordinate  transformation.  The  operator  can  supply  any  mirror 
curvature  to  the  program  as  long  as  all  virtual  image  sources  are 
located  outside  the  resonator.  The  coordinate  transformation  cannot 
treat  internal  focal  points.  It  should  be  pointed  out  that  in  the 
SSCP(ECS),  the  maximum  distance  traversed  in  the  x-direction 
depends  on  the  maximum  value  for  Ax. 

2.  Calculation  Procedure 

The  calculation  procedure  for  the  empty  cavity  configu¬ 
ration  is  as  follows  (see  Fig,  4). 
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Empty  Cavity 


1.  Compute  VR1  and  VR 2  and  VR2  using  eq.  (15). 

2.  Reflect  input  wave  or  waves  from  mirror  1. 
(Maximum  of  six  waves  can  be  specified.  ) 

Mirror  1  and  mirror  2  are  treated  as  if  they 
have  infinite  radius  of  curvature  (flat  mirrors), 
therefore  only  the  spatial  properties  of  the  input 
wave  are  changed. 

3.  Compute  Ax  using  eq.  (17).  (DXMIR1  = 

X WIDTH/NX,  where  XWIDTH  is  an  input 
parameter) . 

4.  Propagate  one  step  to  z  =  Az  using  eq.  (8). 

5.  Compute  Ax  using  eq.  (17). 

6.  Propagate  one  step  in  the  z  direction  using 
eq.  (13). 

7.  Repeat  steps  5  and  6  until  mirror  2  is  reached. 

8.  Reflect  from  mirror  2. 

9.  Repeat  steps  3  through  6  until  mirror  1  is 
reached. 

10.  Store  field  points  on  RSTART  tape  for  use  in 
the  FAR  FIELD  PROGRAM. 

11.  Renormalize  field  amplitude.  This  step  multi¬ 
plies  the  field  at  mirror  1  by  a  constant  so  that 
at  this  location  the  area  under  the  distribution 
remains  fixed  from  pass  to  pass. 

12.  Reduce  the  maximum  spatial  extent  of  the  beam 
from  xmax  to  XWIDTH.  This  step  takes  the  NX 
field  points  which  have  amplitude  and  phase 
information  for  the  beam  spaced  every  DXOUT 
and  replaces  every  field  point  with  beam  infor¬ 
mation  spaced  DXMIR1  apart. 

13.  Repeat  steps  2  through  13  for  NPASS  passes. 

For  the  empty  cavity  configuration,  the  operator  can  determine  the 
total  roundtrip  losses  for  the  resonator  by  summing  the  losses  that 


i 
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occur  at  each  mirror.  The  loss  at  mirror  1  is  determined  by  taking 
the  area  under  the  intensity  distribution  after  completing  step  11  and 
subtracting  from  it  the  area  after  step  2.  The  loss  at  mirror  2  is  the 
difference  in  areas  after  completing  steps  7  and  8  [the  area  after  each 
of  the  indicated  steps  is  printed  as  output  by  both  the  SSCP  and  the 
SSCP(ECS)]. 

The  phase  distribution  obtained  with  the  expanding  coordinate 
system  is  normalized  to  the  geometric  phase  for  the  output  beam.  In 
other  words,  if  there  are  no  perturbations  in  the  resonator,  the  output 
phase  will  be  essentially  flat  (all  geometric  curvature  removed). 

To  include  the  laser  medium  effects  in  the  SSCP(ECS)  an  addi¬ 
tional  step,  other  than  the  steps  outlined  for  the  SSCP,  is  required  to 
compute  the  local  intensities  needed  for  the  gain  subroutine.  Note  in 
Fig.  4  that  propagating  from  mirror  1  to  mirror  2  the  total  extent  of 
the  beam  in  the  x-direction  is  smaller  than  for  the  opposite  traveling 
beam.  Since  the  field  points  of  the  two  opposite  traveling  waves  do  not 
necessarily  coincide,  the  total  local  intensity  (which  is  needed  for  the 
gain  subroutine)  must  be  obtained  by  interpolation.  In  addition  the 
matrix  containing  the  local  intensities  S0T0T  and  GAIN  must  be 
expanded  to  include  field  points  covering  the  entire  beam  of  the  previ¬ 
ous  pass  when  propagating  from  mirror  1  to  mirror  2  at  the  current  z 
location. 

D.  Pseudo  Three-Dimensional  Calculations 

In  high  power  laser  devices  with  transverse  medium  flow,  the 

spatial  homogeneity  of  the  near  field  intensity  distribution  depends 

upon  the  output  coupling  of  the  resonator.  Since  the  coupling  loss  for 

two-  and  three-dimensional  resonators  are  not  equal  [for  two- 

dimensional  confocal  unstable  resonators  the  loss  (1  -  1/M)  is  consid- 

1 

erably  less  than  the  three-dimensional  loss  (1  -  1/M  )],  a  realistic 
near-field  distribution  cannot  be  obtained  with  a  two-dimensional  reso¬ 
nator  calculation.  Furthermore,  three-dimensional  calculations  are 
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necessary  to  determine  output  power  since  this  quantity  also  depends 
upon  coupling  loss.  The  SSCP  and  the  SSCP(ECS)  contain  an  option  to 
perform  a  pseudo-three-dimensional  calculation.  A  description  of  the 
calculation  technique  is  given  below. 

1 .  Calculation  Technique 

The  three-dimensional  cylindrical  resonator  is  divided 
into  two  halves  about  a  line  running  from  the  top  to  the  bottom  of  the 
resonator.  The  same  propagation  algorithms  used  for  the  two- 
dimensional  resonators  can  be  used  if  it  is  assumed  that  for  the  right 
half  (downstream  half)  of  the  resonator  the  field  amplitude  in  the  reso¬ 
nator  is 


=  y3d(r)ej4>  =  “2d 


(*) 


(18) 


for 


7T 

2 


<  6 


0  < 


<  DIA2 
2 


and  for  the  left  half  (upstream  half)  of  the  resonator  for  field  amplitude 
is 


y3d(r’*> =  y3d(r)eJ*  =  y3d(x) 


(19) 


for 


TT 
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and 


-PI  A  2 
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<  x  <  0 
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Thus,  for  the  upstream  half  of  the  resonator,  the  mode  distribution  is 
independent  of  <j>,  and  a  similar  assumption  is  made  about  the  down¬ 
stream  portion.  Next  the  flux  density  that  the  gain  medium  senses 
must  be  adjusted  in  order  to  account  for  the  differences  in  coupling 
loss.  To  find  this  adjustment  factor,  note  that  the  power  contained  in  a 
thin  strip  Ax  wide  and  DIA2  high  in  a  two-dimensional  resonator  is 

Ap^  =  *  DIA2  ■  Ax  ,  (20) 

i  2 

where  I  =  Uj  The  power  contained  in  a  thin  strip  Ar  wide  and  cov¬ 
ering  half  a  three-dimensional  resonator  (either  the  upstream  or 
downstream  half)  is 

Ap^  ,  =  aSR  •  I- ,  •  rAr  =  ASR  •  I,,  •  xAx 
1  3d  3d  3d 

ASR  =  tt  (cylindrical  resonator) 

ui) 

ASR  =  4  (square  resonator) 

ASR  =  0  (two-dimensional  calculation) 

where  the  center  of  the  resonator  mirror  is  assumed  to  be  located  at 
x  =  0,  and  I  is  the  flux  density  in  the  laser  cavity.  Because  of  the 
imposed  angular  requirements  (cylindrical  resonator)  on  in  the 
cavity,  the  three-dimensional  cavity  flux  density  can  be  set  equal  to 
the  two-dimensional  cavity  flux  density.  Thus 

A»3d  =  5D&  ■  x  •  Ap2d  <22> 

However,  in  order  to  satisfy  eq.  (22)  the  flux  density  that  the  medium 
senses  1^  (1^  is  the  flux  density  used  to  calculate  the  gain  which  in  this 
case  is  different  from  the  cavity  flux  I)  must  be  adjusted  as  follows 


(23) 


p  _  ASR 
Lid  "  DIA2  ' 


2.  Calculation  Procedure 

The  procedure  for  computing  the  pseudo-three- 
dimensional  near-field  in  the  presence  of  a  laser  medium  is  similar  to 
the  two-dimensional  procedure  described  for  the  SSCP  and  the 
SSCP(ECS).  The  three-dimensional  procedure  requires  two  additional 
steps.  First,  since  eqs.  (7)  and  (13)  conserve  the  two-dimensional 
area  the  field  amplitude  at  each  propagation  step  must  be  renormalized 
in  order  to  conserve  the  three-dimensional  area.  Each  field  point  is 
multiplied  by  the  following  constant: 


RATIO  =  Renormalization  Constant  (3-d) 


Z|U.m|2|*.  -  x  |Ax 
_i _ 

S|U.m+1|2|x.  -  x  | Ax 
.1  '  l  c 

i 


1/2 


(24) 


where  x  is  the  x-dimension  for  the  center  of  the  resonator.  The 
c 

accumulative  ratio  RCUM3D  is  the  NZ  products  of  RATIO  for  each 
propagation  step  (R/vTIOj  •  RATIO^  RATIO^z) .  Second,  the  total 
local  intensities  must  be  adjusted  according  to  eq.  (23)  before  comput¬ 
ing  the  local  gain.  It  should  be  pointed  out  that  the  local  intensities  are 
not  permanently  changed  by  this  adjustment;  they  retain  their  original 
values  once  the  local  gain  has  been  computed. 

Both  programs,  the  SSCP  and  the  SSCP(ECS)  are  set  up  so  that 
when  ASR  =  0,  two-dimensional  calculations  are  made. 


E.  Program  Uses 

The  Step-by-Step  Computer  Program  (SSCP)  and  the  Step-by- 
Step  Computer  Program  with  Expanding  Coordinate  System 
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[SSCP(ECS)] ,  when  used  to  analyze  empty  cavity  resonators,  will  solve 
for  the  amplitude  and  phase  of  the  lowest  loss  transverse  mode,  the 
total  roundtrip  loss  per  pass,  the  output  coupling  fraction,  the  ampli¬ 
tude  and  phase  of  the  far-field  distribution,  and  the  far-field  normal¬ 
ized  area  distribution.  This  latter  distribution  gives  the  total  area,  or 
power,  under  the  far-field  distribution  for  a  given  far-field  angle. 

Both  near-  and  far-field  information  is  provided  regardless  of  the  res¬ 
onator  configuration.  The  resonator  may  be  stable  or  unstable,  per¬ 
fectly  aligned  or  misaligned  as  a  result  of  mirror  tilt,  or  have  partially 
transmitting  mirrors. 

The  empty  cavity  analysis  is  important  because  it  gives  the 
program  user  the  nonperturbed  intensity  cistributions  for  any  two- 
dimensional  resonator  configuration.  It  also  permits  the  user  to 
determine  the  effects  of  mirror  misalignment  and  coupling  on  the 
intensity  distributions. 

With  an  active  medium  present,  additional  information  can  be 
obtained  from  the  SSCP  and  SSCP(ECS).  Using  the  uniformly  saturable 
gain  routine  contained  in  the  two  programs,  the  user  can  determine 
how  the  intensity  distributions  are  perturbed  when  saturable  gain  is 
present.  Furthermore,  he  can  determine  if  a  given  small  signal  gain 
will  cause  multimoding  in  the  resonator.  For  gain  routines  such  as  the 
Lockheed  2-D,  DF-CO^,  mixing  program,  effects  from  flow,  multi¬ 
wavelengths,  and  refractive  index  gradients  can  be  studied.  In  addi¬ 
tion,  of  course,  parametric  studies  involving  concentration, 
temperature,  and  mixing  of  gas  species  can  be  performed  in  order  to 
determine  their  effect  on  the  output  beam. 

F.  Limitations 

The  two-dimensional  Step-by-Step  Computer  Program  (SSCP) 
has  three  basic  limitations.  The  first  limitation  is  that  the  analysis  is 
limited  to  slit  or  two-dimensional  resonators;  this  restrains  the  laser 
medium  from  varying  in  the  y  direction.  Therefore,  lasers  which  have 
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free  expansion  nozzles,  shock  waves,  etc.  ,  cannot  be  accurately 
analyzed  with  a  two-dimensional  analysis.  The  pseudo-three¬ 
dimensional  calculation  procedure  discussed  in  Subsection  D  permits 
the  user  to  eliminate  the  difference  in  coupling  loss  between  the  two- 
and  three-dimensional  resonators;  however,  the  procedure  still 
requires  medium  uniformity  in  the  y  direction. 

The  second  limitation  of  the  SSCP  is  the  sii.e  of  the  Fresnel 
number  of  the  resonator  which  can  be  analyzed.  For  efficient  calcula¬ 
tions  the  number  of  field  points  NX  should  be  less  than  approximately 
300.  From  past  experience  this  requires  Fresnel  numbers  of  60  or 
less.  The  Step-by-Step  Computer  Program  with  Expanding  Coordinate 
System  [SSCP(ECS)]  removes  this  limitation.  However,  in  removing 
the  Fresnel  number  limit  it  adds  a  constraint  that  no  internal  focal 
points  can  exist  inside  the  resonator  (all  virtual  image  sources  must 
b».  outside  the  resonator). 

The  third  limitation  is  the  inability  to  treat  infinitely  sharp 
mirror  edges.  Since  the  mirror  edges  must  be  gradually  tapered,  the 
fine  structure  in  the  near-field  intensity  distribution  is  eliminated. 
However,  the  overall  structure  is  not  degraded  by  the  tapering  of  the 

mirror  edges  and  the  net  result  is  an  averaging  of  the  fine  structure  to 

produce  a  smooth  varying  intensity  distribution  for  the  transverse 
modes, 

G.  Input  Data  Considerations 

Several  input  variables  to  the  program  are  not  set  by  the 
physical  dimensions  of  the  laser  device  being  analyzed.  A  list  of  these 
along  with  a  brief  description  follows: 

Variable  Description 

NX  Number  of  field  points 

NZ  Number  of  propagation  steps  between 

mirrors  * 
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Variable 


Description 


IGSKIP  Number  of  propagation  steps  between 

GAIN  calls 

AMPINT  Amplitude  of  input  wave  (uniform  or  gaus- 

sian  (nTW /L) 


WOMEGA  1/e  radius  of  input  wave  (L) 

OMEGA  1,  2  Truncation  distance  for  mirror  edges  ( L) 
XWIDTH  Maximum  spatial  extent  of  input  wave  (L,) 


NX  must  be  chosen  such  that  there  are  at  least  two  field  sample  points 
within  the  smallest  length  in  the  x  direction  over  which  the  amplitude  or 
phase  of  the  beam  changes.  The  smallest  length  may  be  OMEGA  or  the 
size  of  a  density  variation  such  as  a  shock  wave.  For  single  mode  lasers, 
NX  will  be  determined  primarily  by  the  truncation  distance  of  the  mirror 
edges  (see  discussion  at  end  of  this  section). 

NZ  is  determined  from  the  following  equation: 

K  2 

Az  <  ~  •  Ax  (25) 


where 


Ax 


XWIDTH 

NX 


and 


Az 


ZLEN 

NZ 


Under  most  conditions,  eq.  (25)  is  sufficient  to  determine  Az;  however, 
if  a  laser  medium  is  present  such  that  the  amplitude  or  phase  of  the 
beam  changes  by  more  than  5  or  10%,  then  Az  must  be  reduced  in 
order  not  to  exceed  this  limit. 

IGSKIP  should  be  set  so  that  the  GAIN  subroutine  is  called 
whenever  the  beam  amplitude  changes  by  5  or  10% 

The  amplitude  of  the  input  waves  (uniform  or  gaussian  distribu¬ 
tion)  is  set  by  AMPINT.  For  empty  cavity  conditions  the  field  distribu¬ 
tion  is  normalized  every  pass  in  order  to  make  the  average  intensity  at 
mirror  1  equal  to  1.0;  therefore,  AMPINT  has  no  meaning  in  this 
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case.  For  gain  filled  resonators  the  distribution  is  not  normalized. 
Since  the  final  mode  distribution  must  be  independent  of  the  input  wave, 
the  value  for  AMPINT  is  arbitrary.  For  values  less  than  the  satura¬ 
tion  levels  of  the  medium,  the  SSCP  and  SSCP(ECS)  require  several 

passes  through  the  resonator  before  the  medium  saturates  (gain  equals 

2 

loss).  By  setting  (AMPINT)  approximately  equal  to  the  saturation 
flux  density  of  the  medium,  the  SSCP  and  SSCP(ECS)  require  the  few¬ 
est  number  of  iterations  for  a  solution. 

W OMEGA  is  the  1/e  radius  for  the  input  wave  or  waves.  The 
amplitude  distribution  for  the  input  is  as  follows: 


| U | input  =  AMPINT  exp(-(x  -  X/ZWOMEGA2)  .  (26) 

For  uniform  distributions,  WOMEGA  —  00 . 

The  truncation  distance  for  the  mirror  edges  is  given  by  OMEGA  1 
and  OMEGA2.  The  values  for  OMEGA  should  be  such  that  2  •  OMEGA < 
0.07  •  DIA.  For  example,  if  DIA2  =  10  cm  then  OMEGA  2  S  0.35  cm. 
The  purpose  of  the  upper  limit  on  OMEGA  is  to  insure  that  the  taper 
region  does  not  contribute  to  a  significant  portion  of  the  mirror  sur¬ 
face.  There  is  an  additional  consideration  that  must  be  made  when 
setting  a  value  for  OMEGA.  Significant  truncation  error  terms  in  the 
propagation  algorithm  can  develop  if  there  are  not  sufficient  sample 
field  points  in  the  taper  region;  therefore,  the  operator  must  be  sure 
that  DXMIR  <  OMEGA/ 3.  Decreasing  OMEGA  means  that  DXMIR 
must  be  decreased  proportionately.  However,  DXMIR  can  only  be 
decreased  by  increasing  NX,  which  means  more  computer  storage  and 
computation  time  to  analyze  the  resonator.  Therefore  the  operator 
must  compromise  between  the  sharpness  of  the  mirror  edge  and  com¬ 
putational  efficiency. 

Since  the  step-by-step  propagation  algorithm  is  a  numerical 
difference  method,  the  first  and  last  field  points  in  the  propagation 
algorithm  must  always  contain  zero.  Therefore  in  the  SSCP  and  the 
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SSCP(ECS),  XWIDTH  must  be  large  enough  to  ensure  that  throughout 
the  entire  cavity,  significant  truncation  errors  are  not  introduced  by 
setting  the  end  points  to  zero  (in  both  programs,  the  value  assigned  to 
the  end  points  are  always  zero).  In  the  SSCP  (where  there  is  no 
expansion  of  the  coordinate  system),  XWIDTH  must  be  large  enough  to 
include  both  geometric  and  diffraction  spreading  of  the  beam.  For 
example,  on  the  first  pass  when  the  initial  beam  is  reflected  from 
mirror  1  and  propagated  toward  mirror  2,  XWIDTH  (which  is  equal  to 
NX  ■  Ax)  must  be  larger  than  the  beam  when  it  arrives  at  mirror  2 
(the  beam  size  at  mirror  2  will  depend  on  both  the  radius  of  curvature 
of  mirror  1  and  diffraction  spreading).  In  the  SSCP(ECS),  geometric 
spreading  is  accounted  for  by  the  expanding  coordinate  system. 
Therefore,  XWIDTH  which  sets  the  width  of  the  initial  input  beam  does 
not  have  to  be  made  large  enough  to  include  geometric  spreading,  but 
only  large  enough  to  include  diffraction.  This  is  the  advantage  of  the 
expanding  coordinate  system.  When  using  the  SSCP(ECS),  XWIDTH 
usually  is  slightly  larger  than  DIA1,  while  in  the  SSCP,  XWIDTH 
usually  is  larger  than  the  largest  cavity  mirror. 


i 
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II. 


program  description 


A.  Summary  of  the  SSCP(ECS) 

The  Step-by-Step  Computer  Program  with  Expanding  Coordinate 
System,  or  SSCP(ECS),  is  composed  of  a  variety  of  FORTRAN  sub¬ 
routines  and  functions.  For  simplicity,  they  can  be  grouped  into  four 
divisions  based  on  function:  primary,  satellite,  utility,  and  input/ 
output  routines.  The  boundaries  between  divisions  are  arbitrary;  but 
in  general  reflect  the  importance  or  difficulty  of  implementation  of  a 
subprogram. 

Primary  routines  are  generally  the  longest  and  most  logically 
complex  of  a  system  of  routines.  In  addition,  they  usually  contain  the 
fundamental  mathematical  concepts  and  techniques  upon  which  the  pro¬ 
gram  was  based.  The  SSCP(ECS)  system  consists  of  three  such 
routines . 


PR5 1 

Main  program,  defines  input,  contains 
overall  logic 

P50SBS 

Propagates  one  or  more  waves  through 
free  space 

P50MED 

Propagates  one  or  more  waves  one  way 
through  a  medium 

Satellite 

routines  are  generally  self-contained  routines  which 

accomplish  some  smaller  portion  (often  repetitive)  of  the  overall  solu 
tion.  Four  of  these  routines  are  incorporated  in  the  SSCP(ECS). 

P50GAI 

Returns  a  set  of  gain  values  at  a  particu¬ 
lar  location  in  the  medium 

GAINA 

Computes  the  gain  for  each  element  of 
each  wave  for  some  location  in  the 
medium 

MIRROR 

Reflects  one  or  more  waveforms  from  a 
particular  mirror 

SETMIR 

Sets  up  a  mirror  vector 

29 


The  GAINA  routine  is  described  as  a  satellite  routine  only  because  of 
its  mathematical  simplicity  and  its  short  length.  A  program  such  as 
the  Lockheed  DF-CO^  mixing  program^would  obviously  be  a  primary 
routine  when  substituted  for  GAINA. 

Utility  routines  generally  have  little  to  do  with  the  calculation 
that  the  program  system  is  performing,  but,  accomplish  some  neces¬ 
sary  function  instead.  The  four  routines  of  SSCP(ECS)  in  this  category 
are : 


P50GET 

Retrieves  one  intensity  value  from  the 
last-pass  intensities  matrix  by  a  triple 
linear  interpolation 

P50INT 

Reduces  the  output  wave(s)  from  fully 
expanded  form  to  compacted  form  for 
start  of  next  pass  by  linear  interpolation 

NORM2D 

Two-dimensional  area  conservation 
{renormalization  routine) 

RAT3DA 

Three-dimensional  area  conservation 
routine 

Input/output  routines  are  normally  special  purpose  routines 
designed  to  manipulate  a  particular  type  of  input  or  output  in  a  particu 
lar  way,  and  hence  are  self  explanatory  and  highly  inflexible.  Three 
routines  of  this  class  are: 

AREA21  Computes  and  prints  the  area  under  a 

curve 

PWROUT  Computes  and  prints  the  near  field  power 
and  power  transmitted  for  one  or  more 
waves  positioned  before  the  output  mirror 

RSTART  Restarts  a  previous  run  by  reading  a 

magnetic  tape 

The  fourth  routine  in  this  category,  ALLOUT,  is  in  reality  a 
package  of  routines  designed  to  allow  the  programmer  to  describe  the 
type,  frequency,  and  quantity  of  output  he  desires  while  retaining 
flexibility.  Since  these  routines  have  nothing  to  do  with  the  problems 
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which  the  SSCP(ECS)  was  designed  to  handle,  their  description  will  be 
given  in  Appendix  A. 

13.  Main  Program  and  Key  Subprogram  Descriptions 

1.  PR51  (Main  Program) 

The  main  program,  Pi  ' ,  defines  all  input  parameters 
and  controls  all  output  either  directly  (by  calling  the  output  package)  or 
indirectly  (by  supplying  an  output  file  code  to  routines  which  use  the 
output  package).  In  addition,  PR51  controls  job  flow  by  initializing 
certain  parameters  before  each  case  begins,  defining  certain  vectors 
(such  as  the  mirror  vectors),  and  defining  the  manner  of  execution  of 
the  primary  routines  P50MED  and  P50SBS. 

All  important  arrays  are  declared  in  PR51  and  supplied  to  sub¬ 
programs  through  the  argument  list  or  by  labeled  COMMON.  This 
feature  permits  the  user  to  alter  the  size  of  certain  arrays  as  desired 
to  minimize  the  storage  requirement.  All  arrays  dimensioned  in  PR51 
are  indexed  by  one  or  more  of  five  quantities:  (1)  the  maximum  num¬ 
ber  of  propagating  points  in  the  x-direction,  (2)  the  maximum  number 
of  last-pass  intensity  values  stored  in  the  z-direction,  (3)  the  largest 
number  of  propagating  waves,  (4)  the  largest  number  of  passes  that 
PR51  can  handle,  and  (5)  the  largest  number  of  gain  values  that  PR51 
can  handle.  These  five  values  are  assigned  to  variables  NXDIM, 
NZDIM,  NWDIM,  NPDIM,  and  NGDIM,  respectively,  and  appear  in  a 
DATA  statement  at  sequence  number  1820.  Table  I  describes  the  rela¬ 
tionship  between  each  subscript  of  each  array  and  the  approximate 
quantity  involved.  The  abbreviations  used  in  the  table  are 

X  s  NXDIM 
Z  s  NZDIM 
W  =  NWDIM 
P  =  NPDIM 
G  -  NGDIM 
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TABLE  1.  PR51  ARRAY  SUBSCRIPT  REFERENCES 


Array  Name 

Array  Type 

T 

Complex 

U 

Complex 

V 

Complex 

XMIR1 

Complex 

XMIR2 

Complex 

GAIN 

Complex 

SO 

Real 

SONEW 

Real 

SOTOT 

Real 

XLAMDA 

Real 

XKZERO 

Real 

AMPINT 

Real 

INFOUT 

Real 

WOMEGA 

Real 

X 

Real 

OUT1 

Real 

OUT  2 

Real 

OUT3 

Real 

OUT4 

Real 

OUT5 

Real 

OUT6 

Real 

OUT7 

Real 

Subscript  Reference 
No. 


Should  the  user  need  to  change  one  or  more  of  the  dimensions  of  PR51, 
he  should  first  identify  which  arrays  are  affected  (by  Table  I)  and  then 
change  each  array  declaration  in  PR51  to  reflect  the  changed  dimen¬ 
sion.  Last,  he  should  change  the  appropriate  value  in  the  DATA 
statement  at  sequence  number  1820  so  that  all  subprograms  will  auto¬ 
matically  handle  the  change. 

2 .  Selected  Subprogram  Descriptions 

Flowcharts  and  listings  are  provided  in  Appendix  E  for 
each  primar\  routine  and  the  more  complex  satellite  and  utility  rou¬ 
tines.  Other  routines,  including  the  output  package  (GPOP),  have  only 
listings.  Documentation  describing  the  purpose  of  each  routine  and  the 
nature  of  the  arguments  supplied  are  provided  in  the  listings  them¬ 
selves  with  the  exception  of  three  GPOP  routines.  The  description  of 
GPOP  usage  appears  in  Appendix  A.  The  information  provided  by  the 
listings  (and,  if  present,  the  flowchart)  is  generally  sufficient  to  ade¬ 
quately  describe  the  subprogram.  However,  three  routines  warrant 
further  mention. 

P50GAI 

The  critical  feature  of  P50GAI  is  the  storage  and 
retrieval  of  the  last-pass  intensities  prior  to  the  call  to  the  gain  rou¬ 
tine.  Since  a  major  portion  of  dimensioned  storage  lies  in  the  SO 
(last-pass  intensity)  matrix,  a  facility  exists  by  which  the  user  can 
avoid  storing  intensities  at  every  medium  propagation  step.  When 
retrieving  intensity  values  from  the  matrix,  therefore,  a  linear 
interpolation  procedure  is  used  to  derive  an  approximate  intensity 
value  for  each  field  point.  As  the  propagation  proceeds  along  the 
z-direction,  current  intensity  values  are  stored  in  place  of  last-pass 
intensities  which  are  no  longer  needed  so  that  the  wave  re-entering  the 
medium  during  the  next  P50GAI  call  will  contain  new  values.  The  SO 
index  in  the  z-direction  increments  when  propagating  toward  the  back 
mirror,  and  decrements  when  propagating  toward  the  output  mirror. 
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SETMIR 


The  mirror  setup  routine,  SETMIR,  computes  three 
components  for  each  location  of  the  mirror  matrix.  The  first, 
reflectivity,  is  supplied  directly  to  the  routine  as  an  input  parameter 
but  may  be  locally  adjusted  by  SETMIR  to  accommodate  edge  tapering 
(see  Section  I.  B.4).  The  second  component,  mirror  curvature,  is 
defined  as 


CURV.  =  exp[-jKo(x.  -  xc)2/RCURV]  , 


i  =  1,  2,  •  •  •  ,  NX 


where 

x  is  the  x-value  at  the  mirror  center 

c 

RCURV  is  the  mirror  radius  of  curvature. 
The  third  component,  mirror  tilt,  is  defined  as 


j[2K  p(x  -x  )] 

TILT.  =  e  °  1  c  ,  i  =  1,  2,  •  •  •  ,  NX 


where 


0  =  mirror  tilt  angle  (radians). 

GAINA 

The  gain  routine  returns  an  array  defining  the  power 
gain  and  refractive  index  variation  of  the  medium  for  each  wave  and 
each  field  point  at  a  desired  z-location.  The  gain  matrix  is  complex, 
with  the  real  part  containing  the  local  power  gain,  gj,  and  the  imagi¬ 
nary  part  containing  the  refractive  index  variation,  Ar|.KQ.  The 
GAINA  routine  supplied  with  the  SSCP  and  SSCP(ECS)  only  returns  a 
power  gain. 
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c. 


Input  Description 


PR51  is  the  only  routine  in  the  SSCP(ECS)  package  requiring 
input  to  be  prepared  by  the  user.  All  such  input  is  entered  on  data 
ca  -ds  by  the  FORTRAN  NAMELIST  convention,  with  which  the  reader 
is  assumed  to  be  familiar.  (A  complete  description  can  be  found  in  the 
FORTRAN  reference  manual.)  In  general,  a  NAMELIST  data  record 
consists  of  a  "$"  in  column  2  of  the  data  card  followed  immediately  by 
the  particular  NAMELIST  name  for  this  record,  followed  by  one  or 
more  blanks  or  a  comma,  followed  by  the  list  of  variables  and  corre¬ 
sponding  values  (in  the  form  of  assignment  statements  separated  by 
commas),  and  finally  terminated  by  a 

Two  NAMELIST  group  names  are  used  by  PR51,  the  first, 
CASENO,  used  to  define  the  number  of  cases  to  be  analyzed  during  the 
run,  and  the  remaining  one  used  for  each  case  to  describe  all  run 
conditions,  geometry,  and  constraints.  The  data  deck  for  a  run  will 
always  have  as  its  first  NAMELIST  record  a  card  of  the  form 

$CASENO  NCASES  =  5$ 

to  describe  the  number  of  cases  to  be  analyzed  (here  set  to  "5"  for 
illustration).  Following  this  card  will  be  NCASES  groups  of  cards, 
each  group  having  a  first  card  starting  "$INPUT.  .  .  "  and  a  last  card 
ending  with  a  "$"  following  the  last  data  value  defined. 

The  parameters  which  must  be  defined  in  the  case  description 
block  named  "INPUT"  are: 

(Title) 

CASE  Real  number  in  format  XXX.  XX  to  be 

displayed  for  titling  purposes  on  output 
as  case  number 
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(Run  option  descriptions) 


NPASS 

IFILE 

IRCODE 

IRTAG 

IROUT 


The  number  of  passes  which  the  starting 
wave  is  to  be  propagated  (no  upper  limit) 

The  FORTRAN  file  code  on  which  all 
normal  output  will  appear  (normally  =  6) 

FORTRAN  file  code  for  restarting  a 
previous  run,  if  >0  (must  match  the  file 
code  assigned  to  a  magnetic  tape  on  a 
system  control  card) 

Pass  number  from  which  restart  informa¬ 
tion  is  desired  from  file  code  IRCODE,  if 
IRCODE  >  0 

FORTRAN  file  code  for  generating  a  tape 
for  later  restart  or  far  field  analysis,  if  >0 
(must  match  the  file  code  assigned  to  a  mag¬ 
netic  tape  on  a  system  control  card) 


(Cavity  description) 

NX  Number  of  propagating  field  points  (see 

Section  I-G  and  II-B-1) 

XWIDTH  Minimum  spatial  extent  in  x  of  the  beam 

(see  Fig.  4),  L* 

NZ  Number  of  propagation  steps  to  be  taken 

between  mirrors  (see  Section  I-G  and 
II-B-1) 

ZJLEN  Distance  between  mirrors,  L 


XCENTR  Relative  location  along  x  of  mirror  axis 

centerline,  L 


ASR 


Pseudo -three -dimensional  parameter 
(see  Section  I-D-l) 


#Unit  of  length  is  arbitrary  but  must  be  consistent  for  a  data  case. 
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(Medium  description) 

ZMDLOC  Location  of  the  center  of  the  medium 
relative  to  the  output  mirror,  L  (see 
Fig.  2) 

ZLENMD  Length  of  the  medium,  L  (see  Fig.  2) 

(Gain  routine  description) 

IGSKIP  Skip  parameter  for  call  to  GAINA,  where 

IGSKIP  =  3  would  call  the  gain  routine 
every  third  propagation  step  (see  Sec¬ 
tion  I-G) 

GZERO  gQ,  the  small  signal  gain  used  by 

GAINA,  L"1 

XIZERO  I0,  the  saturation  flux  density  used  by 

GAINA,  watts/L^ 

(Input  wave  description) 

NWAVES  The  number  of  propagating  waves 

XLAMDA  \,  or  wavelength  for  each  wave;  array  of 

NWAVES,  L 

AMPINT  Initial  amplitude  for  each  wave  (see  Sec¬ 

tion  I-G);  array  of  NWAVES,  (watts)  1/2/L 

WOMEGA  w  for  each  gaussian  input  wave,  if  >  0  (see 
Section  I-G);  array  of  NWAVES,  L 

(Mirror  descriptions) 

DIA1  Diameter  of  output  mirror,  L 

DIA2  Diameter  of  back  mirror,  L 

REFL1  Power  reflectivity  of  output  mirror* 

REFL2  Power  reflectivity  of  back  mirror* 

RCURV1  Radius  of  curvature  of  output  mirror,  L 

RCURV2  Radius  of  curvature  of  back  mirror,  L 

*Fraction  with  values  between  0  and  1. 
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OMEGA  1  1/e  truncation  distance  for  output  mirror 

edges  (see  Section  I-G),  L 

OMEGA2  1/e  truncation  distance  for  back  mirror 

edges  (see  Section  I-G),  L 

TILT1  Tilt  angle  of  output  mirror,  radians 

TILT2  Tilt  angle  of  back  mirror,  radians 

TLOSS1  Power  loss  of  output  mirror* 

(Output  option  selection) 


INFOUT  Near  field  print  options  (a  "zero"  stored 

in  the  ith  location  of  this  vector  will  sup¬ 
press  the  output  for  the  matching  pass 
number;  a  "one"  will  produce  the  normal 
output);  array  of  NPASS 


D.  Restart  Facility 

PR51  controls  the  generation  of  a  restart  tape  if  a  certain 
option  is  supplied,  while  RSTART  controls  the  search  of  a  restart  tape 
to  find  and  read  a  desired  waveform  if  two  other  parameters  are 
supplied.  The  user  will  normally  not  be  concerned  with  the  format  of 
the  data  on  the  tape  since  the  program  handles  the  information  solely 
on  the  basis  of  three  input  variables  during  case  definition  (IRCODE, 
IRTAG,  and  IROUT).  Should  it  be  necessary  to  know  the  format  of  the 
data  on  the  restart  tape,  the  following  information  should  be  sufficient. 

The  information  on  the  restart  tape  is  written  by  FORTRAN 
unformatted  (binary)  WRITE  statements.  The  first  record  consists  of 
five  words  regarded  as  titling  information  only.  The  second  record 
consists  of  20  words  used  for  describing  the  quantity  of  information  to 
be  read  from  the  tape.  The  first  word  of  the  integer  array  is  NX,  the 
second  is  NZ,  and  the  third  is  IWAVE,  where  these  variables  repre¬ 
sent  the  number  of  propagating  field  points  in  the  x-direction,  the 
number  of  SO  matrix  locations  used  for  storage  in  the  z-dimension, 


^Fraction  with  values  between  0  and  1. 
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and  the  number  of  propagating  waves  (NWAVES  in  data  input  section), 
respectively. 

The  remaining  logical  records  on  the  tape  occur  in  groups  of 
three,  each  group  representing  the  results  of  a  different  round-trip 
pass. 

Record  Information  Contained 

1  The  pass  number  for  this  set  of  infor¬ 

mation.  (This  number  is  searched  for  if 
a  restart  is  requested.) 

Z  The  near-field  pattern  for  each  propa¬ 

gating  wave,  written  and  read  by 
( [U(I,  J) ,  1=  1 ,  NX] ,  J=  1 ,  IWAV E) 

3  The  last-pass  intensities  matrix,  so 

that  propagation  can  be  continued; 
written  and  read  as 
[([SO(I,  J,K),  1=1,  NX],  J=l,  NZ), 

K=  1,  IWAVE)] 

Groups  of  three  logical  records  are  read  at  a  time  in  the  event  of  a 
restart  until  the  pass  number  read  from  the  tape  matches  the  restart 
number  supplied  by  input  variable  IRTAG. 

E.  Gain  Routine  Replacement 

GAINA,  the  gain  routine  provided  with  SSCP(ECS)  is  designed 
to  be  replaced  with  a  more  sophisticated  routine  whenever  desired.  A 
replacement  of  this  nature  will  be  simplest  if  the  new  routine  uses  the 
same  name  and  calling  arguments  as  does  the  present  GAINA  routine. 
If  additional  arguments  (such  as  gQ  and  Iq  in  GAINA)  are  required, 
provision  for  their  initialization  can  be  made  by  including  a  labeled 
COMMON  block  in  both  PR51  and  the  gain  subroutine  and  initializing 
the  parameters  by  a  READ  statement,  an  initialization  subroutine,  or 
by  a  BLOCK  DATA  subprogram.  Modification  of  the  argument  list 
required  for  the  gain  routine  is  not  recommended  since  at  least  the 
P50GAI  routine  would  have  to  be  changed,  and  possibly  P50MED,  to 
assure  that  the  arguments  are  passed  successfully  to  the  gain  routine. 
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III.  MACHINE  UTILIZATION 

The  SSCP(ECS)  package  is  written  in  FORTRAN  IV.  The 
program  was  originally  written  for  a  GE635  computer  operating  under 
the  GECOS  III  operating  system,  though  this  manual  reflects  changes 
made  allowing  the  program  to  run  on  a  CDC  6600  computer  operating 
with  the  SCOPE  monitor  (Version  3.4).  This  requires  about  34K 
(decimal)  or  103K  (octal)  words  of  memory,  of  which 

NXDIM[(NZDIM  +  7)NWDIM  +  12]  +  NWDIM(3  •  NGDIM  +  6) 

+  NPDIM  +  NZDIM 

resides  in  dimensioned  arrays  in  the  main  program.  For  the  current 
dimensioning  parameters 

NXDIM  =  101 

NZDIM  =  21 

NWDIM  =  2 

NPDIM  =  101 

NGDIM  =  1001 

the  total  storage  required  for  dimensioned  variable  storage  is 
13,000  words  (decimal)  or  31,300  words  (octal). 

Processor  time  depends  on  too  many  parameters  and  option 
settings  to  be  accurately  estimated,  but  the  longest  case  analyzed  so 
far  consumed  about  two  minutes  on  the  CDC  6600  and  about  15  minutes 
on  the  GE635. 
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APPENDIX  A 


GENERAL  PURPOSE  OUTPUT  PACKAGE  (GPOP)  DESCRIPTION 

The  GPOP  package  was  written  to  eliminate  the  need  for 
specialized  output  routines  (and  to  eliminate  the  inevitable  errors 
incurred  during  their  development)  while  still  retaining  enough  flexi¬ 
bility  that  output  could  be  tailored  to  a  program's  needs.  The  routine 
is  able  to  handle  real  or  complex  vectors  and  perform  various  normal¬ 
ization  operations  on  the  data  to  be  displayed.  Tables  and/or  printer 
plots  can  be  specified  for  any  vector  as  can  a  skip  parameter  in  case 
the  data  need  not  be  tabled  in  their  entirety.  Pages  are  titled  and  data 
are  headed  by  the  appropriate  user-supplied  column  heading.  Up  to 
seven  vectors  of  information  can  be  handled  at  once.  Since  all  such 
output  is  fully  specified  by  providing  the  appropriate  information  in  a 
FORTRAN  subroutine  CALL  statement,  the' possibility  of  user  error  is 
virtually  eliminated. 

Usage 


CALL  ALLOUT  (IFILE,  NPTS,  ISKIP,  NMAT,  NCHAR,  LABEL, 
ITYPEj.IOPTj.NCHARj,  LBL^AMATj,  DUMMY  r 


ITYPE  ,  IOPT  ,  NCHAR  ,  LBL  ,  AMAT  .DUMMY  ) 
n  n  n  n  n  n 


where  n  =  NMAT  and  1  <  n  <  7. 
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Description  of  Parameters 


for 
each 
vector' 


IFILE 

Output  file  code 

NPTS 

Number  of  points  in  each  user  vector 

ISKIP 

Skip  parameter  for  vectors  (i.e.  ,  use 
every  ISKIPth  location) 

NMAT 

Number  of  vectors  supplied  to  ALLOUT,  <  7 

NCHAR 

Number  of  characters  in  title 

LABEL 

Title  (either  an  array  or  a  Hollerith 
literal) 

ITYPE. 

1 

Matrix  type  (1  =  real,  2  =  complex)  for 
ith  vector 

JOPT. 

i 

th 

Output  option  (described  below)  for  i 
vector 

NCHAR. 

l 

Number  of  characters  in  label  for  i*h 
vector 

LI3L. 

i 

th 

Label  for  i  vector  (array  or  Hollerith 
literal) 

AMAT. 

i 

Name  of  user  matrix  to  be  displayed 

DUMMY. 

l 

Dummy  vector  supplied  by  user  (dimen¬ 
sioned  at  least  NPTS) 

The  IOPT^  option  is  a  seven-digit  integer  which  completely 
describes  the  sort  of  output  to  be  produced  from  the  corresponding 
vector.  Each  digit  has  a  specific  meaning.  For  the  description  below, 
assume  that  the  integer's  digits  are  numbered  from  right  to  left  (i.  e.  , 
7654321),  and  that  numbers  shorter  than  seven  digits  will  have  zeroes 
assumed  for  the  leftmost  digits. 
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Digit  Number 


Options 


Z 


3 


Value  type  to  be  taken  from 
vector 

0  number  is  assumed 

real 

1  real  part  of  complex 

number 

Z  imaginary  part  of 

complex  number 

3  magnitude  of  complex 
number 

4  intensity  of  complex 
number 

5  pha.se  (radians)  of 
complex  number 

6  phase  (degrees)  of 
complex  number 

Option  for  tabular  output 

0  no  tabular  output  for 

this  vector 

1  produce  tabular  output 

for  this  vector 

Option  for  printer  plot 

0  no  plot  for  this  vector 

1  plot  vector  versus 
AMATj,  full  page 
scaling 

2  plot  vector  versus 
AMATj,  plot  scaled  for 
"convenient"  axis 
increments 
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Options 

Not  currently  used 
Option  for  normalization  type 
0  no  normalization 

1  vector  normalized  to 
one  by  division 

2  vector  normalized  to 
zero  by  subtraction 

6  Option  for  position  of  normalization 

0  no  normalization 

1  center  point  is  normal¬ 
ized  to  zero  or  one 

2  highest  value  point  nor¬ 
malized  to  zero  or  one 

3  lowest  value  point  nor¬ 
malized  to  zero  or  one 

7  Option  for  normalized  area  plot 

0  produce  vector  itself 

1  use  vector  to  produce 

normalized  area  plot 


Digit  Number 

4 

5 
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2785-8 


ALLOUT 


ONEOUT 


VALUE 


LBLPOS 


main  routine  of  the  GPOP  set;  produces  a 
complete  set  of  listings  and/or  plots  when 
referenced  by  a  CALL  statement 

prepares  one  vector  to  be  listed  and/or 
plotted;  called  once  for  every  vector 
supplied  to  ALLOUT 

computes  the  value  for  one  location  in  the 
vector  to  be  listed  and/or  plotted;  called 
for  each  location  of  every  vector  handled 
by  ONEOUT 

positions  alphanumeric  column  heading 
above  each  vector  to  be  listed 


PGPLT2  produces  line  printer  plot 

CSCALE  scales  axes  for  PGPLT2 


Fig.  A1 .  GPOP  subprogram  references. 
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APPENDIX  B 


FAR  FIELD  PROGRAM  DESCRIPTION 

I.  SYSTEM  DESCRIPTION 

A.  Introduction 

The  purpose  of  the  FAR  FIELD  PROGRAM  (FFP)  is  to 
determine  the  far-field  intensity  distribution  of  the  laser  output  beam, 
the  total  far-field  power,  and  the  incoherent  average  of  the  near-  and 
far-field  intensity  distributions  (averaging  is  only  necessary  if  the 
laser  is  operating  multimode). 

The  method  used  to  determine  the  far-field  intensity  distribu¬ 
tion  is  to  take  the  near-field  distribution  that  is  transmitted  by  the 

output  mirror  and  multiply  each  field  point  by  a  phase  correction  fac- 

2 

tor  exp  (i  Kq(x  -  x^)  /(RCURV),  where  xc  is  located  at  mirror  center. 
This  correction  factor  acts  as  if  the  output  beam  is  reflected  from  a 
spherical  mirror  having  a  radius  of  curvature  RCURV .  This  is  referred 
to  as  a  far-field  mirror.  The  far  field  is  by  definition  located  in  the 
focal  plane  of  a  spherical  mirror,  therefore  propagating  the  beam 
reflected  from  the  far-field  mirror  a  distance  ZLENFF  =  RCURV/2, 
the  far-field  intensity  distribution  is  obtained  for  the  output  laser 
beam . 

Incoherent  averaging  is  used  whenever  the  output  beam  is 
multimode  or,  in  other  words,  whenever  the  transverse  mode  does  not 
stabilize  but  continues  to  change  from  pass  to  pass.  The  averaging  is 
performed  by  adding  the  intensities,  at  each  field  point,  to  each  suc¬ 
cessive  pass  through  the  resonator  and  dividing  by  the  number  of 
passes  averaged.  For  multimode  lasers,  both  the  averaged  near-  and 
far-field  intensity  distributions  will  converge  to  a  stable  solution  after 
sufficient  averaging. 
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B.  Calculation  Procedure  (Empty  Cavity) 

For  two-dimensional  (ASR  =  0)  far-field  distributions  the 
computational  procedure  is  straightforward.  The  near  field  amplitude 
and  phase  for  the  empty  cavity  is  obtained  from  the  RSTART  tape. 

This  distribution  is  multiplied  by  mirror  1  (output  mirror)  and  the 
resulting  distribution  (e.  g.  ,  the  portion  of  the  beam  that  spills  past  the 
edges  of  the  mirror  plus  the  portion  transmitted  by  mirror  1)  multiplied 
by  the  far-field  mirror.  Finally  the  beam  is  propagated  to  the  focal 
plane  of  the  mirror  by  a  two-dimensional  Fresnel  integral  (see  Ref.  1 
for  a  description  of  the  Fresnel  Integral).  The  user  must  choose 
RCURV  so  that  the  spatial  extent  of  the  far  field  distribution  is  similar 
to  the  spatial  extent  of  the  near  field.  For  example,  if  RCURV  is  too 
small,  the  number  of  sample  points  containing  significant  information 
will  be  a  small  percentage  of  the  total  number  of  sample  points  NX. 

When  ASR  t  0  or  in  other  words  when  pseudo  three-dimensional 
calculations  are  being  made,  the  far-field  calculation  procedure  is 
similar  to  that  described  above  except  the  two-dimensional  Fresnel 
integral  must  be  replaced  with  a  three-dimensional  Fresnel  integral. 
The  following  propagation  algorithm  is  used  to  find  the  far-field  distri¬ 
bution  (in  a  transverse  plane  cutting  through  the  focal  point)  for  a 
cylindrical  near-field  distribution  in  Section  I-D. 


U(x) 


iK 


Nx/2 


"r2 

j=l 


iK 

e 


(x2+{jAx)2) 

>  RCURV- 


•  U(jAX)AX  - 


Nx 

y 


iK 

e 


(xZ+(jAx)2) 
o  RCURV 


j=Nx/2 


(  /  KQxjAx  v  /  KqxJAx  V\ 

\ o\ZLENFf  /  1  o\ZLENFF  JJ 


U(jAx)  Ax  , 


(1) 


for  x  <  0  , 
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where  x  is  measured  in  the  far  field  plane  with  the  center  located  at 
x  =  0,  U(jAx)  is  the  output  beam  (includes  both  real  and  imaginary 
terms),  and  is  the  Struve  function.  For  x  >  0,  the  far  field  distri¬ 
bution  is  found  by  changing  the  sign  in  front  of  the  Struve  functions. 

C .  Calculation  Procedure  (Laser  Medium  Included) 

Since  the  step-by-step  propagation  algorithm  requires  that 
mirror  edges  be  tapered,  the  phase  distribution  for  the  near  field 
transverse  mode  will  change  significantly  for  small  changes  in  x  in  the 
region  of  the  taper  (note  that  this  rapid  phase  change  occurs  in  the 
taper  region).  This  change  is  seen  as  an  increase  in  the  phase  curva¬ 
ture  as  | x J  increases.  In  empty  cavity  resonators,  the  amplitude  of 
the  beam  in  the  taper  region  falls  rapidly  and  therefore  the  small 
amplitude  combined  with  a  rapidly  changing  phase  does  not  significantly 
distort  the  far-field  distribution.  However,  in  transverse  flow  lasers 
where  the  upstream  portion  of  the  resonator  can  have  considerably 
higher  unsaturated  gain  (especially  in  the  taper  region),  the  amplitude 
of  the  beam  in  the  taper  region  may  no  longer  be  insignificant  and 
therefore  there  could  be  significant  distortion  in  the  far  field  due  to  the 
taper  region.  To  correct  for  this  effect  which  is  introduced  because  of 
the  necessity  to  taper  the  mirror  edges,  there  is  an  option  in  the  FAR 
FIELD  PROGRAM  to  take  the  phase  of  the  output  beam  from  the  gain 
filled  resonator  and  subtract  from  it  the  phase  of  the  empty  cavity 
output  beam  (obtained  from  empty  cavity  RSTART  tape).  This  removes 
the  phase  curvature  due  to  the  taper  region,  but  does  not  affect  any 
phase  distortion  introduced  by  the  medium. 

The  use  of  the  phase  subtraction  is  an  operator  option.  Experi¬ 
ence  has  shown  that  the  phase  curvature  in  the  taper  region  for  the 
SSCP(ECS)  is  considerably  smaller  than  for  the  SSCP  and  therefore 
phase  subtraction  has  not  been  needed  when  the  SSCP(ECS)  was  used. 

To  better  understand  the  phase  subtraction,  refer  to  Fig.  Bl.  Fig¬ 
ure  Bl(a)  shows  a  typical  near-field  phase  distribution  for  an  empty 
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Fig.  Bl.  Example  of  phase  subtraction 
when  computing  far  field 
intensity  distributions  in  the 
far  field  program. 


confocal  unstable  resonator.  The  phase  is  essentially  flat  in  the  center 
of  the  resonator  and  increases  rapidly  in  the  taper  region.  The  phase 
distribution  for  the  medium  filled  resonator  looks  essentially  the  same 
as  the  empty  cavity  distribution  except  for  some  distortion  near  the 
center.  The  difference  in  the  two  distributions  is  shown  in  Fig.  Bl(c) 
and  is  essentially  flat  everywhere  except  for  the  distortion  caused  by 
the  medium.  The  far  field  distribution  for  the  last  two  phase  distribu¬ 
tions  (assuming  some  amplitude  distribution)  will  be  the  same  if  the 
amplitude  is  small  in  the  taper  region.  However,  if  the  amplitude  is 
not  small  in  the  taper  region,  the  far-field  distribution  for  the  second 
phase  distribution  will  appear  distorted.  Thus  by  using  the  third  phase 
distribution,  the  rapidly  varying  phase  in  the  taper  region  is 
eliminated . 
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II. 


PROGRAM  DESCRIPTION 


A.  Summary  of  the  FFP 

The  Far  Field  Program  (FFP)  is  composed  of  primary, 
satellite,  and  input/ output  routines  of  the  sort  described  earlier.  The 
two  primary  routines  are 

FFLD  Main  program,  defines  input,  contains 

overall  logic 

FFLDAV  Propagates  one  or  more  waves  to  the  far 

field;  controls  averaging;  controls  output 
of  far  field  and  averaging  information 

Four  satellite  routines  are  contained  in  the  FFP. 

PMAT30  Sets  up  a  Fresnel  integral  vector  for 

two-dimensional  propagation  to  the  far 
field 

SETM01  Sets  up  the  far  field  mirror  vector 

FFPRP2  Propagates  a  wave  to  the  far  field  by 

three-dimensional  propagation 

STRUVE  Computes  the  Struve  function  by  Taylor 

series  approximation 

The  three  input/output  routines,  RSTART,  AREA21,  and 
ALLOUT  have  been  discussed  in  Section  II- A  of  the  SSCP(ECS)  part  of 
this  manual  and  will  not  be  further  discussed  here. 

Four  of  the  routines  deserve  further  mention  at  this  point. 

FFLD 

FFLD,  the  main  program,  controls  the  definition  of  the  input 
near  field  wave  prior  to  calling  the  far  field/averaging  routine 
(FFLDAV).  The  procedure  in  constructing  the  near  field  wave  is  to 
define  both  a  "primary"  and  a  "secondary"  wave,  where  the  latter  may 
not  be  supplied  for  certain  cases.  When  phase  subtraction  of  the 
empty  cavity  waveform  is  desired,  the  secondary  wave  is  restored 
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from  the  specified  magnetic  tape  and  the  phase  correction  terms 
computed.  As  each  primary  wave  is  restored  (from  a  different  tape), 
it  is  phase-corrected  by  the  secondary  wave  (if  it  has  been  supplied)  to 
yield  the  desired  near  field  waveform. 

There  is  no  functional  difference  between  restart  tapes  gener¬ 
ated  by  the  SSCP  and  those  generated  by  the  SSCP(ECS)  —  both  are 
fully  compatible  with  the  tape  requirements  of  the  FFP  package.  The 
difference  lias,  in  the  past,  been  that  an  SSCP  tape  will  also  require 
the  empty  cavity  secondary  tape  for  phase  subtraction  while  the 
SSCP(ECS)  primary  tape  exhibits  a  sufficiently  flat  phase  profile  to  not 
warrant  the  empty  cavity  tape.  The  latter  situation  is  a  significant 
advantage  since  the  number  of  tapes  should  be  minimized  for  faster  job 
turnaround. 

SETM01 

The  mirror  setup  routine,  SETM01,  computes  three  components 
for  each  location  of  the  mirror  matrix.  The  first  is  the  transmission 
distribution.  The  distribution  takes  into  account  the  size  of  the  output 
mirror  DIA  1  ,  the  output  mirror  reflectivity  REFL  which  is  constant  over 
the  mirror  surface,  and  the  transmission  loss  TLOSS1  of  the  output 
mirror.  This  distribution  when  multiplied  by  the  near  field  distribution 
taken  from  the  RSTART  tape,  gives  the  field  distribution  incident  on  the 
far-field  mirror.  The  second  component,  curvature,  is  defined  as 

CURVi  =  exp  [-jKo(x.  -  xc)2/RCURV],  i  =  1,2,  •••,  NX 

The  third  component,  mirror  tilt,  is  defined  as 

j[2K  6(x  -x  )] 

TILT.  =  e  °  1  c  ,  i  =  1,2,  •  •  •  ,  NX 

where 

p  =  mirror  tilt  angle 

x  =  x  at  the  mirror  center 
c 

RCURV  =  radius  of  curvature  of  the  far-field  mirror 
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PMAT30 


PMAT30  generates  the  Fresnel  Integral  matrix  by  which 
diffraction  propagation  is  performed.  To  simplify  computation,  the 
matrix  is  symmetrical  around  the  (2-NX-l)  point.  In  general,  the  ith 
value  of  the  matrix  is  defined  as 


fiZT  jKQ[(i-I)Ax]2/2L 
Ax  •  v/zLX  e 

where 

L  s  ZLENFF,  or  far-field  propagation  distance 

though  the  procedure  in  PMAT30  performs  this  computation  in  a  dif¬ 
ferent  form  than  that  shown. 

FFLDAV 

The  concept  of  averaging  is  included  in  the  FFLDAV  routine 
along  with  the  far-field  propagation  (which  is  straightforward). 
Averaging  can  be  used  independently  on  either  the  near-field  or  the 
far-field  intensity  distribution  or  both,  and  consists  of  a  numerical 
average  at  each  point  of  all  the  values  of  the  intensity  distribution  over 
some  set  of  passes. 

B.  Input  Description 

FFLD  is  the  only  routine  in  the  FFP  package  requiring  input  to 
be  prepared  by  the  user.  All  such  input  is  entered  on  data  cards  by  the 
FORTRAN  NAMELIST  convention,  with  which  the  reader  is  assumed  to 
be  familiar.  (A  complete  description  can  be  found  in  the  FORTRAN 
Reference  Manual.  )  In  general,  a  NAMELIST  data  record  consists  of 
a  "$"  in  column  2  of  the  data  card  followed  immediately  by  the  particu¬ 
lar  NAMELIST  name  for  this  record,  followed  by  one  or  more  blanks 
or  a  comma,  followed  by  the  list  of  variables  and  corresponding  values 
(in  the  form  of  assignment  statements  separated  by  commas),  and 
finally  terminated  by  a  "$.  " 
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Two  NAMELIST  group  names  are  used  by  FFLD,  the  first, 
CASENO,  used  to  define  the  number  of  cases  to  be  analyzed  during  the 
run,  and  the  remaining  one  used  for  each  case  to  describe  all  run  con¬ 
ditions,  geometry,  and  constraints.  The  data  deck  for  a  run  will  always 
have  as  its  first  NAMELIST  record  a  card  of  the  form 

$CASENO  NCASES  =  5$ 

to  describe  the  number  of  cases  to  be  analyzed  (here  set  to  "5"  for 
illustration) . 

Following  this  card  will  be  NCASES  groups  of  cards,  each 
group  having  a  first  card  starting  "$INITL.  .  .  "  and  a  last  card  ending 
with  a  "$"  following  the  last  data  value  defined. 

The  parameters  which  must  be  defined  in  the  case  description 
block  named  "INITL"  are: 


A  SR 

Pseudo -three -dimensional  parameter 
(see  Section  I-B  in  this  appendix) 

IPSTRT 

Pass  number  at  which  to  initially  restore 
the  primary  wave  (>0) 

IPFIN 

Last  desired  pass  for  which  a  near  field 
wave  is  to  be  constructed  (sNPASS  of  run 
generating  tape) 

ISTAG 

The  pass  number  of  the  wave  to  be  desig¬ 
nated  the  secondary  wave  (sNPASS  of  run 
generating  tape) 

IPCODE 

Primary  wave's  file  code,  >0 

ISCODE 

Secondary  wave's  file  code,  if  >0 

IFILE 

The  FORTRAN  file  code  on  which  all 
normal  output  will  appear  (normally  =  6), 
>0 

N  WAVES 

The  number  of  waves  to  be  taken  from 
tape  (<NWDIM  in  FFLD) 

ZLENFF 

The  far  field  propagation  distance 
(normally  half  the  far  field  radius  of 
curvature) 

RCURV 

The  far  field  radius  of  curvature  (see 

Section  I-B  in  this  appendix) 

XKO 

The  wave  number  (2tt/\.)  for  each  wave; 
array  of  NWAVES,  L“1 

DX 

The  spacing  between  field  points 

NX 

The  number  of  field  points  (iNXDIM  in 
FFLD) 

TLOSS 

The  power  loss  of  the  output  mirror* 

REFL 

The  power  reflectivity  of  the  output 
mirror* 

BETA 

Tilt  angle  of  far  field  mirror,  radians 

DIA1 

Diameter  of  output  mirror,  L** 

DIAZ 

Diameter  of  back  mirror,  L 

VR2 

Virtual  source  distance  for  the  back 
mirror;  (set  to  1.E30  if  SSCP(ECS)  not  used) 

ZLEN 

Distance  between  cavity  mirrors,  L 

INFOUT 

Near  field  print  options;  array  of  IPFIN 
(<101) 

th 

0  stored  in  the  i  location  of  this 

vector  will  suppress  the  near  field 
output  for  the  matching  primary  pass 
number 

1  stored  in  the  i^  location  of  this  vec¬ 

tor  will  produce  full  near  field  output 
for  the  matching  primary  pass 
number 

^Decimal  percent  between  0  and  1 
**Units  are  arbitrary  but  must  be  consistent  for  a  case. 
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IFFOPT 


IFFOUT 


Far  field  propagation  and  near/far  field 
averaging  switch,  where  the  i^h  location 
in  IFFOPT  controls  the  matching  primary 
pass  number 

0  Do  not  propagate  to  FF;  do  not 
average  NF  or  FF 

1  Propagate  to  FF;  do  not  average  NF 
or  FF 

2  Do  not  propagate  to  FF;  average  NF 
only 

3  Propagate  to  FF;  average  NF  only 

4  Propagate  to  FF;  average  FF  only 

5  Propagate  to  FF ;  average  both  NF 
and  FF 

Far  field  averaging  output  options,  where 
the  i**1  location  in  IFFOUT  controls  the 
output  for  the  matching  primary  pass 
number 

0  Print  neither  set  of  results 


1  Print  FF  results  only 

2  Print  averaging  results  only 

3  Print  both  FF  and  averaging 
results 
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C .  Flowcharts  and  Li  stings 
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Sample  Case  for  FFP 


The  case  described  by  the  input  below  is  a  far-field  analysis  of 
the  waveforms  computed  during  the  sample  case  "8.45"  described  in 
Appendix  D.  A  magnetic  tape  was  generated  and  saved  by  that  program 
and  will  be  used  for  input  to  the  FFP.  The  user  is  assumed  to  be 
familiar  with  the  conventions  adopted  by  his  own  computer  center  when 
referencing  a  previously  saved  magnetic  tape. 

The  far-field  analysis  is  to  be  two-dimensional  (ASR  =  0.) 
adhering  to  the  same  mode  of  analysis  used  in  the  SSCP(ECS)  sample 
case  8.45.  The  tape  is  to  be  read  from  its  beginning,  rather  than 
some  intermediate  pass  number,  and  is  to  read  no  further  than  the 
7^  pass.  One  wave  is  to  be  taken  from  the  tape  at  each  pass  (the  tape 
could  contain  more  waves).  The  file  code  for  the  primary  wave  is 
logical  file  code  10  (the  user  is  assumed  to  have  previously  informed 
the  computer  system  that  file  code  10  is  to  be  associated  with  the 
SSCP  restart  tape).  There  is  no  secondary,  or  empty  cavity,  tape  to 
furnish  phase  terms  for  subtraction  from  the  primary  waveforms.  The 
file  code  for  normal  program  output  is  logical  file  code  6.  The  far- 

4 

field  propagation  distance  is  to  be  6  x  10  cm.  The  wave  number 
(XKO)  can  be  found  in  the  fifth  NAMELIST  group  (WAVES)  at  the  begin¬ 
ning  of  the  normal  output  of  the  SSCP(ECS)  for  case  8.45  following  the 
vector  name  XKZERO.  The  tlx  mesh  spacing  can  similarly  be  found  in 
the  second  NAMELIST  block  (CAVITY)  following  the  term  DXOUT. 

The  number  of  mesh  points  along  the  x  direction  is  known  (from  the 
input  to  8.45)  to  be  81.  The  first  seven  passes  read  from  the  restart 
tape  are  to  have  their  near  field  values  printed,  to  be  propagated  to  the 
far  field,  and  to  have  the  far  field  values  printed.  No  averaging  is  to 
be  performed.  The  output  mirror  of  the  8.45  cavity  is  100%  reflecting 
(the  loss  has  no  meaning  in  this  instance  since  no  light  energy  passes 

through  the  mirror).  The  radius  of  curvature  of  the  far  field  mirror  is 

4 

twice  the  propagation  length  above,  or  12  x  10  cm.  The  far-field 
mirror  is  not  tilted.  The  diameters  of  the  output  mirror  and  back 
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mixTor  for  the  8,45  cavity  are  5.74  cm  and  11.  1  cm,  respectively. 
The  virtual  radius  of  the  back  mirror  is  found  in  the  second,  or 
CAVITY,  NAMELIST  group  of  case  8.45  after  the  term  VR2.  The 
distance  between  mirrors  in  the  cavity  is  10,000  cm. 

Sample  Input  Deck 


(column  2) 


$CASENO  NCASES=  1$ 

$INITL 

ASR=0.  ,  IPSTRT=0,  IPFIN=7,  NWAVES=1, 

IPCODE=10,  ISCODE=0,  ISTAG=0, 

IFILE=6,  ZL,ENFF=6.  E4,  XK0(  1)  =  20944.  , 

DX  =  .  225,  NX  =  8 1 , 

INFOUT(l)  =  7*l,IFFOPT(l)=7*l,IFFOUT(l)=7*l, 
TLOSS=0.  ,  REFL=  1 .  ,  RCURV=1. 2E5,  BETA=0.  , 

DIA  1  =  5. 74,  DIA2=  1 1 .  1,  VR2=1.  E30, 

ZLEN=  10000.  $ 

The  output  produced  by  the  FFP  consists  of  four  major  groups. 
The  first  is  not  optional  and  includes  the  following: 

1 .  Repetition  of  the  input  data 

2.  Description  of  the  far-field  mirror  matrix 
(intensity  and  phase  listings  and  plots) 

3.  Restart  summary  of  titling  information  con¬ 
tained  on  the  secondary  tape,  followed  by  a 
listing  and  plot  of  the  phase  correction  informa¬ 
tion  (appears  only  if  a  secondary  wave,  or  empty 
cavity,  tape  is  used) 

4.  Description  of  the  "ideal"  uniformly  illuminated 
aperture  for  this  cavity  configuration 

5.  Far  field  summary  of  the  wave  in  4. 

6.  Restart  summary  of  titling  information  contained 
on  the  primary  tape 

The  near -field  summary  (optional)  lists  and  plots  the  intensity 
and  phase  information  of  a  wave  taken  from  a  primary  tape  (after 
phase  correction,  if  a  secondary  tape  was  specified). 


The  far-field  summary  (optional)  lists  and  plots  magnitude, 
intensity,  normalized  magnitude,  normalized  intensity,  normalized 
phase,  and  normalized  area  of  a  near-field  wave. 

The  averaging  summary  (optional)  lists  and  plots  the  incoherent 
average  of  both  near-field  and  far-field  waveforms. 

The  following  pages  represent  the  output  produced  by  the  sam¬ 
ple  data  case.  All  types  of  output  except  the  secondary  wave  restart 
and  averaging  are  represented. 
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TIOSS  0,  REFL  0 ,  1C0D000  0E  01  RCURV  0:11000000E  06  RET A  0. 

DI*1  0,S74ooflOOE  01  DU2  O.UIOOOOOE  02  VR?  O-IBOOOOOOE  3l  ZLEN  0.10000000E  05 

END  N*M6UST  INITL 


MIRROR  niScRlRTtON 


I  x 

1  ^0 . 90000000E 

2  «0',87750e00i 

3  eO i 8550SOOOE 

4  s0i83249499E 

5  =0 . B1000000E 

6  -0 , 787900001 

7  -0, 7ft5f>0800E 

8  =0 . 7429O0DOE 

9  sO . 770008006 

10  i0.69750800E 

11  =0 . 67500000E 

12  =0,612S0800E 

13  -0.63000BOOE 

1*  s0.6fl750000E 

15  s0.505o0800E 

16  g0,56250B00E 

17  =Q,54000000E 

18  =0.5l750000E 

19  =0,49500800E 

20  s0.47250000E 

21  =  0 . 49000000E 

22  ct  42750000E 

23  s0.405n0800E 

24  iO  3S250B00E 

25  e0.38000000E 

26  i0.33750800E 

27  •0,315o0800E 

26  «0.29250600E 

29  c0,270fi0800E 

30  «0,24750600E 

31  cO,22SeoeooE 

32  aO;2fl250000E 

33  iO.lBOOOOOOE 

34  .-.15750800E 

35  «0.13500800E 

36  c0.1l250800E 

37  iO.OnOOOOOOE 

38  sO  67000000E 

39  *0  45000800E 

40  =0,2?500800E 

41  0, 

42  0  22500600E 

43  0.45000000E 

44  0.675000001 

45  Ql9nOOOOOOE 

46  0.11250000E 

47  0.13500000E 

48  0.15790000E 

49  O.lOOQOOOOE 

50  0.20290000E 


INTERS 

01 

O.IQQDOOOOE 

01 

0 i 1Q000000E 

01 

0  •  lOOQOQOOE 

01 

0.99Q99996E 

01 

0.99999996E 

01 

o.iooonoooE 

01 

0 i lOOQOOOOE 

01 

o.iooonoooE 

01 

0.10000000E 

01 

OilOOOOOOOE 

01 

0 .lOOQOOOOE 

01 

0.99999998E 

01 

0 . 99999998E 

01 

0 . 99999999E 

01 

0 .lOOQOOOOE 

01 

0.99999T98E 

01 

0.99999999E 

01 

0.99999999E 

01 

0.10000000E 

01 

o.iooaooooE 

01 

0.99999999E 

01 

0 . 10000000E 

01 

0.99999999E 

01 

0  '■  99999999E 

01 

0 . 10030000E 

01 

O.IOOOOOOOE 

01 

0.10000000E 

01 

0. 

01 

0. 

01 

0. 

01 

0. 

01 

0. 

01 

0, 

01 

0. 

01 

0, 

01 

'1, 

00 

0. 

00 

0. 

00 

0, 

00 

0. 

0. 

00 

0, 

00 

0, 

00 

0. 

00 

0. 

01 

0, 

01 

0, 

ot 

0. 

01 

0. 

01 

0. 

I  Tv  PHASE  (DO) 

01  -0  *.  90002237F.  0? 
01  •0.30008240E  02 
01  *0) 11026748E  02 
fin  0 J 26942243E  02 
00  0.63898710E  02 

01  0199842685E  02 

01  0J13477415E  03 

ni  0jl68693llE  03 
01  *0 *  18840182E  03 
01  «0 j 12650767E  03 
01  =0  95626431E  02 
00  io  65757497E  02 
00  sO . 36901069E  02 
00  a0'.90571454E  01 
01  0 J 17774270E  02 

00  0.43593178E  02 

00  0 J 60399585E  02 

00  0  9Z.93480E  02 

01  0  11497487E  03 

01  0 J 1367*376E  03 

UO  0 i 1I750013E  03 
01  0:i77244giE  03 

00  «0|l6402fi02E  03 
00  -0;i4630715E  03 
01  >0  129608B0F  03 
01  »0 ■ 11390899E  03 
01  »0 J 99225614E  02 
iO  50000800E  00 
-0;50000800E  00 
iO  50000000E  00 
iO,5DOOOBOOE  00 
iO)500S060DE  DO 
•0 , 50000000E  00 
iO.SOOOOOOOE  00 
iO  .500008 OOE  00 
iOjSOOOOOOOE  00 
sO  50000000E  00 
io'.sooooeooe  oe 
iO  15OOOO0QOE  00 
iO.SOOOOOOOE  00 
.0  50000800E  00 
ioisoooooooE  00 
iO'.SOOOOOOOE  00 
iolsoooooooE  00 
iOjSOOOOOOOE  30 
io  50000800E  00 
iO.SOOOOOOOE  00 
iojsooooeooE  oo 
>0)50000000E  00 
iO iSOOOOOOOE  00 
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APPENDIX  C 


SSCP  DESCRIPTION  AND  USAGE 

I.  PROGRAM  DESCRIPTION 

A.  Summary  of  the  SSCP 

The  Step-by-Step  Computer  Program  (SSCP)  shares  a  number 
of  attributes  with  the  SSCP(ECS)  program  described  earlier  in  this 
manual.  As  an  illustration,  three  of  four  satellite  routines,  two  of  four 
utility  routines,  and  all  four  input/output  routines  of  the  SSCP(ECS) 
program  are  used  by  both  programs.  The  three  routines  which  are  not 
shared  are  peculiar  to  the  SSCP(ECS)  package  and  bear  no  equivalents 
in  the  SSCP  solution.  All  three  of  the  primary  routines  are  different, 
as  one  might  suspect  from  the  definition  of  a  primary  routine,  though 
their  logic  structure  is  strikingly  similar. 

The  three  primary  routines  used  by  the  SSCP  are 

PR23  Main  program,  defines  input,  contains 

overall  logic 

PRPSBS  Propagates  one  or  more  waves  through 

free  space 

PRPGAI  Propagates  one  or  more  waves  one  way 

through  a  medium 

The  three  satellite  routines  shared  with  SSCP  (ECS)  are 

GAINA  Returns  a  set  of  gain  values  at  a  particular 

location  in  the  medium 

MIRROR  Reflects  one  or  more  waveforms  from  a 

particular  mirror 

SETMIR  Sets  up  a  mirror  vector 
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The  two  utility  routines  shared  with  SSCP(ECS)  are 


NORM2D 

Two-dimensional  area  conservation 
(renormalization)  routine 

RAT3DA 

Three-dimensional  area  conservation 
routine 

‘  \ 

The  input/output  routines  are  exactly  the  same  as  in  th<6 
SSCP(ECS) ,  namely 

AREA21 

Computes  and  prints  the  area  under  a 
curve 

PWROUT 

Computes  and  prints  the  near  field  power 
and  power  transmitted  for  one  or  more 
waves  positioned  before  the  output  mirror 

RS  TART 

Restarts  a  previous  run  by  reading  a 
magnetic  tape 

ALLOUT 

The  general  purpose  master  subroutine 
for  the  GPOP  routines 

B.  Description  of  Primary  Routines 

1.  PR23  (Main  Program) 

All  important  arrays  are  declared  in  PR23  and  supplied 
to  subprograms  through  the  argument  list  or  by  labeled  COMMON. 

This  feature  permits  the  user  to  alter  the  size  of  certain  arrays  as 
desired  to  minimize  the  storage  requirement.  All  arrays  dimensioned 
in  PR23  are  related  to  one  or  more  of  four  quantities:  (1)  the  maxi¬ 
mum  number  of  propagating  points  in  the  x-direction,  (2)  the  maximum 
number  of  last-pass  intensity  values  stored  in  the  z-direction,  (3)  the 
largest  number  of  propagation  waves,  and  (4)  the  largest  number  of 
passes  that  PR23  can  handle.  These  four  values  are  assigned  to 
variables  NXDIM,  NZDIM,  NWDIM,  and  NPDIM,  respectively,  and 
appear  in  a  DATA  statement  at  sequence  number  1860.  Table  2 
describes  the  relationship  between  each  subscript  of  each  array  and 
the  appropriate  quantity  involved.  The  abbreviations  used  in  the  table 
are 


1 16 


X 


NXDIM 


Z  -  NZDIM 
W  s  NWDIM 
P  =  NPDIM 

Should  the  user  need  to  change  one  or  more  of  the  dimensions  of  PR23, 
he  should  first  identify  which  arrays  are  affected  (by  Table  II)  and  then 
c'  i.ge  each  array  declaration  in  PR23  to  reflect  the  changed  dimen¬ 
sion,  Last,  he  should  change  the  appropriate  value  in  the  DATA 
statement  at  sequence  number  I860  so  that  all  subprograms  will  auto¬ 
matically  handle  the  change. 

2.  Primary  Subprogram  Descriptions 

Only  two  subprograms  are  not  shared  with  the  SSCP(ECS), 
these  being  the  primary  subroutines  responsible  for  propagation 
through  free  space  and  propagation  through  a  medium.  Their  logic 
structure,  however,  is  almost  indistinguishable  from  theirECS  coun¬ 
terparts,  differing  mainly  in  the  different  equations  used  in  propaga¬ 
tion.  For  this  reason,  no  further  description  of  these  routines  will 
follow  other  than  a  listing  for  each. 

C.  Input  Description 

The  SSCP  and  SSCP(ECS)  share  exactly  the  same  input  method, 
with  all  but  two  variables,  NX  and  XWIDTH,  retaining  their  values 
independent  of  the  program  selected.  The  considerations  to  be  made 
when  determining  values  for  these  two  variables  are  presented  in 
Section  I-G  of  the  SSCP(ECS)  User's  Manual. 
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TI. 


MACHINE  UTILIZATION 


The  SSCP  package  is  written  in  FORTRAN  IV.  The  program 
was  originally  written  for  a  GE635  computer  operating  under  the 
GECOS  III  operating  system,  though  this  manual  reflects  changes 
made  allowing  the  program  to  run  on  a  CDC  6000  computer  operating 
with  the  SCOPE  monitor  (Version  3.4).  This  requires  about  41k 
(decimal)  or  1Z7K  (octal)  words  of  memory,  of  which 

NXDIM[(NZDIM+8)NWDIM+13]+NPDIM+5NWDIM 

resides  in  dimensioned  arrays  in  the  main  program.  For  the  current 
dimensioning  parameters 


NXDIM  =  301 

NZDIM  =  21 

NWDIM  =  2 

NPDIM  =  101 

the  total  storage  required  for  dimensioned  variable  storage  is 
21, 480  words  (decimal)  or  52,  000  words  (octal). 

Processor  time  depends  on  too  many  parameters  and  option 
settings  to  be  accurately  estimated. 
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III. 


LISTINGS 
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APPENDIX  D 


SAMPLE  CASE  FOR  SSCP(ECS) 


The  case  (named  "8.  45")  described  by  the  input  data  below 
involves  10  round-trip  passes  with  restart  information  to  be  written  on 
file  code  10  (a  control  card  must  appear  elsewhere  in  the  deck  assign¬ 
ing  file  code  10  to  a  new,  or  scratch,  magnetic  tape)  and  normal  pro¬ 
gram  output  to  be  directed  to  file  code  6.  The  input,  or  starting,  wave¬ 
form  for  this  case  is  not  to  be  taken  from  a  restart  tape  but  is,  instead, 
to  be  generated  by  parameters  described  below.  The  number  of  mesh 
points  along  the  x-axis  (across  the  face  of  a  mirror)  is  to  be  81,  and 
the  minimum  width  over  which  these  81  points  are  to  be  defined  is 
9  cm  (see  Fig.  4).  There  are  to  be  101  mesh  points  (100  steps)  between 
mirrors,  where  the  mirrors  are  10,  000  cm  apart.  The  x-value  of  the 
mirror  centerline  is  to  be  0.  The  analysis  will  be  two-dimensional 
(ASR  =  0.).  The  medium,  whose  center  is  located  5000  cm  from  the 
output  mirror,  is  to  be  5000  cm  in  width,  thereby  filling  half  the 
cavity.  Tne  gain  routine,  to  be  called  at  every  propagation  step  in  the 
Z -direction,  is  supplied  two  parameters  ~  gQ  and  IQ>  with  values 
3.  5  x  10  **  and  100.  ,  respectively.  One  wave  of  wavelength  3  x  10  ^  cm 
and  uniform  initial  amplitude  1.0  is  to  be  propagated.  The  initial  wave 

does  not  have  a  gaussian  structure.  The  output  mirror  has  a  diameter 

4 

of  5.  74  cm,  100%  reflectivity,  a  radius  of  curvature  of  -2  x  10  cm, 

an  edge  taper  of  0.45  cm,  no  tilt,  and  no  loss  in  the  mirror  itself 

(meaningless  in  this  case  since  the  reflectivity  was  total).  The  back 

mirror  has  a  diameter  of  11.1  cm,  100%  reflectivity,  a  radius  of 

4 

curvature  of  4  x  10  cm,  an  edge  taper  of  0.  6  cm,  and  no  tilt.  Normal 
output  is  desired  for  all  10  passes. 


Sample  Input  Deck 


(column  2) 

$CASENO  NCASES=  1$ 

$INPUT 

CASE-8.  45,  NPASS=  10,  IFILE  =  6,  IRCODE=0,  IROUT=10, 
IRTAG=0,NX  =  81,XWIDTH  =  9.  ,NZ  =  101,  ZLEN=10000.  , 

XCENTR  =  0.  ,  ASR=0.  ,  ZMDLOC  =  5000.  ,  ZLENMD  =  5000.  , 
IGSKIP=  1 ,  GZERO=3 .  5E-4,  XIZERO=100.  ,  NWAVES=1, 
XLAMDA(1)=3.  E-4,  AMPINT(1)=1.  ,  WOMEGA(  1)=0.  , 

DIA  1  =  5.  74,  REFLU1.0,  RCURVl  =  -2.  E4,  OMEGA  1  =  .  45, 

TILT  1  =0.  ,  T LOSS  1=0.  ,  DIA2=11.  1 ,  REFL2=  1.0,  RCURV2=4.  E4, 
OMEGA2=.  6,  TILT2  =  0.  ,  INFOUT(  1  )=  1 0*1  $ 

The  output  produced  by  either  the  SSCP  or  SSCP(ECS)  for  the 
data  above  consists  of  three  major  groups.  The  two  that  are  not 
optional  are  the  summation  of  the  data  input  defining  the  case  and  the 
summary  information  printed  during  each  round-trip  propagation.  The 
first  includes  the  following  information: 

1.  Repetition  of  the  information  read 

2.  Descriptions  of  both  mirror  matrices  (intensity 
and  phase  listings  and  plots  for  each) 

3.  Description  of  each  input  wave  (magnitude, 
intensity,  phase,  and  normalized  area  listings 
and  plots  for  each). 

The  round-trip  summary  information  is  self-explanatory, 
basically  monitoring  the  area  of  the  propagating  wave  as  it  proceeds 
through  a  round-trip  pass. 

The  third  type  of  output  (optional)  that  can  be  produced  summa¬ 
rizes  the  waveform  itself  at  the  conclusion  of  a  round-trip  pass.  The 
magnitude,  intensity,  normalized  magnitude,  normalized  intensity, 
normalized  phase,  and  normalized  area  are  listed  and  plotted  versus 
their  corresponding  x-location. 

The  following  pages  represent  the  output  produced  by  the  sam¬ 
ple  data  case.  All  three  types  of  output  are  represented. 
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INPLT  W*VF  DtSrB;pTtON 
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TO  3000006*0  TO  TO  3SS»S>S»Z*0  00  3Z9T9T9T9*0-  TO  3T60606Q»*0-  TO  3*9£9£9Cl’0» 

TO  3»9£9C9£**0  TC  3T60U060»*0  00  32«T9T9T9‘0  TO  SSSfrSfrifrZ'O"  TC  ilZilLcLf T8  300000006 


starting  p«ss  numpe»  i  for  case  a. as 

AREA ( S I  AT  START  qF  PASS 

2D  AREA  0. 1?775006  02 

after  interpolation 

2D  AREA  0.9U2S90E  01 

AFTER  OUTPUT  MIRROR 

2D  AREA  0.5T389V1E  01 

AFTER  ST£p»BY=STfcP  PROPAGATION 
2D  AREA  0.573*6096  01 

AFTER  MEOIUM 

2D  AREA  0.3206506E  02 

AFTER  STEP. Br-STEP  PROPAGATION 
2D  AREA  0.3207932E  02 

AFTER  BACK  MIRROR 

20  AREA  0<3O14538E  02 

AFTER  STEP-BY-STEP  PROPAGATION 
2D  AREA  0 , 3014541E  02 

ATTER  MEDIUM 

2D  AREA  g;iSl6766E  03 

AFTER  ST6P5BY-STEP  PROPAGATION 
2D  AREA  0.1916770E  03 

FINISHED  PASS  NUMBER  1 

POWER 

near  field  power*  0.1516770E  0.» 

OUTPUT  PUHEP  APTfcR  MIRROR*  0.6655180E  02 

POWER  OlFPEPkNCE*  0.8512518E  02 
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NE*R  r 1 E L n  P0W£»*  0.1777500E  02 

OUTPUT  POWgR  *FT£R  H1RR0R*  0,67500006  01 

POWER  DIFFERENCE*  0 . 1102500E  02 
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Results  rap  qquui  thip 


HOC  □  O  O  o  o  O  O  OO  OOOCOOOOOiaoOOCOOOOOaOOOOOOHrtHHOoOOQOOO 

ooo  o  oooooo  O  CD  oo  a  000000000000000  aoooooooooooooooo  ooo 

III! 

UUUJU^JUUUUULJJ  UiJLiUUUuluJLJULJuiUUuiUUiULiUiUiUlJUUUuiULJULiltJaiUJLJLUJU 
c  a:  n  (v  h  ho  4  cvi  our-  p*.  c.  4  to  am  hnn  ^  do  «©  <e  o  <r  rv  jn  h  n  hCj  r)  ^  r)  a  .h  <©  4  k-  cv  k-  .h  <r>  4 

«t  DIMS  Ht  t  in  NO  HMHlDHHHOOf  HHB)  a  c  4  ro  OO  Ifl  (O  O  WO  CV  C  O  OHniOO  o  ,hc  K>0  C\J  O'  o 

>.,0  NHCO^ff  onrO«JjMO(Mr  .H  C  HUM^X.  NCHflN^HHOttOHH-T.(Vfl  H  0(\J 
a  oo  c  o  c  r  m  c  co  o>  o-  <oo  it  chi  <  fo  if,  <3  ru  t  ru  h  r  nco  <0  c  bj  cr  in  ho  o  hit  r  a  c  o  o  c  <c 

<  oooooo-or-  4  k-  id  4  ru  co  o-o  hno  on  nin  t  h  n  r  o  op <0  c  o  m  o  «  or*  a  c  k*  oBioinooON 

00  c*  c  o-  o  ?  o  o  ®  r  o  ^  aiN«c  of«o»N(vi3aHHNNhx  Hfo  <r  ;ic  t  Nncoorom  t  a  rtt^  10  ho 

ooo'ooo  c  c  c  c  c  cr  oo®r’OintMOlfiH1oH<JHOH<o{M®inHrfo®*  009  9  oot  <o  no  hia  co 

VI  HOOOO  >  C  C  O  O  C  C  OOOOO-OOffltOfflrr-CCJX^^  ^KJfOPOPJf^HH  H<3  HHO  T-#r4»H(MfM»OK>PO 

Q  •.•••»••■•  ••.•  •  ».  •  ,  •«  *.  •  .  ..  «,  •.  »  •  •.  a.  •  «. 

o  c  o  o  c  o  o  o  o  o  a  ao  ocaococcocooocooocooococoocoocooococoo 


OHHH^HrtfiHHHHOOOHOOOOOOQrtHJMHHHHHHNHrtftlHHHrt 
ooocoooooaoo  ooooooooocr  00000c  000000000000 
*  •  it  it  (i  n  »  i»  *  1  u  i»  o  n  11  11  <i  n  11 

ui  u.1  uj  ai  iu  ui  a.'  xuiu:  uj  uj  u.'  lu  uj  Ui  ui  u  u:  uj  uj  u;  m  4J  to  ui  m  hi  ^juj  l^i  gi  u;  iuuj  UJ  m  >iKtMii 

iii  I  I  (\  O  C  t>  H®  f  (P  in  <  fOfVlAr’OK)f')fOMC«^^<C'CO,CO'ro  ri  ^  C  o  H  O  HIT  ri 

1 ft  ®  fo  00  *)  ao  r  00  o>  «m  Ki  oi>  t  ^  oj\r<o«K)iAAjN«innoro®(\jp)OHHMm 

<1  in  at  nr  o>  o  o  ho  n  cn  c:  <;  k  aj  cv  00  cNr  th^poitxopoolkj  qsj  it  »  r  4  o  a;  c  in  r* 

X  r-irioo^oD'Cr-  ofvfo®  o*b  00  oc  ♦  hit  »tocmoo®ocPs®NP)oc®®o 
a  r-  o  MBtirniMoaoP)  4  ooro  no 

n  r  r  (\j4flr  foa)  in  r*  00  00  ooeio®NrtONoto>toHo  00  n  o  >o  t  o>  r  in  n 

•  r  -OK  K5  HQCOCK!NHin®NM'JHfPOO  N»«HCO®0-»  H  ^  K  <C  |A  (Vt  H  «H4>  f\| 

:.  v  h  cm  cj  .h-o  humcj  *4r4r*r>»-taoc.ir>rototvtv  ho  cv r>>  ro  r*  c  c>  k  pjcj  h  *4  cmm  in  to  h 

ar  . . . . . 

o  c^oaocociaacD  ocoocor,c  iOocoo3coocoocaocoooOoocOooao 

Z  #  I  I  lit  •••••••••••  l|  •• 

■C  m  4  4  r  f\j  (V  JVj  H  H  H  r-HHOOOCDODO^CCOCOC  COCCGCOOCOCOCOCCC  oo 
ooo  OO  OOOOOO  OOOOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOOO 

u.u;u:uju:uju;u-'u;ljiu  uju.uju.u<a^:u.ujujuju;u-u-ujuju;ujujuiu.ujujuju;ujujujujaju.ujujuj  u-  uj  u:  lu 

ck  a  n  o  cvr*  r-  4  rvm  hiA  4  H«*oro«»  two  ok  ao  ro  4  fMo*  k  «  hooia  o  onro  h  h«  (ni  (vn  t  k 

•  (MIA  4  OOO  ONONOHf  TOKJp  CV  O  C  C  O^OVHK®«(MlT  HOO>|n®«ClAHOH^ClA(MA®0 

b-  ®  ®  f  iA  o  o  4  c  r*  04  «ino  hpo«  K)  oc*  4  not  .h  cv  m  k  4  m  .h  t>  4  ro  t  ano«««  on  o  4  n  t  o 

mm  «-t  tv  c  4  cv  m  Hin  <o  ooco  <vm  do  hhuao  ro  -o  hh<m>  h  ^  cvk  co  rg  or  for  orato  r  cm  4  h 

—  (\j  c  o-opocucnojo  rv<.  H««NO*Hflr  10  oe  o»«  ««  no*«  o«ir.  ®  nonjcrinr  cwcro® 

IT.  C  HO>  oONBJMf.ONtlA  OOC  t  CM®  «  (\JO  HIAK  OK  P)  tH  C  H  4  CV  4  O  O  CO  M  HM  Q  «  □  4  (Nl  4  H 

•  r)  ^  4  H  H  HH(®  HKK  K3tA»lANHt»tOKOMA04<OatHKJKNTt(MO«KCON4,4WKW 

X  ♦  N  HK  fO  HW  -O  H  *-l(M  4  4)0  HWKJfAK  «  O  H»»»»COKKKKKgO(QflOlCCOKKKOCO(0(OQKK 

X  . . . . . 

o  oooooooooooo  oooooooooooooooooocoooooooooooooooooo 

Z 


H  H  H  H  f\i  H  H<M  H 

ocoo  o  a  ooo 

<1  (1  It  II  II  II  II  II  (I 

LtJ  til  Ui  UJ  UJ  UJ  UJ  UJ  U.’ 

oooo  ®  cv  r  ire 

04  h<o  too  4  m  in 

«>  0>  C  K  Ht  K  <CIA 

o®  4i®  ®  fo  tv  m 
o  nr)  o  c  4  a  in 
poiako*  4  ro  oro 

NO  HH  MT  ®  K 
HMlrttt  CJH  HN 


r 

a 


rXMNIMHHHHCOC  OOCOOOOO 

ooooooooooo  ooo  oo  ooo 

I  I  I  »  •  •  I  • 

uj  ui  aJ  ul  uj  uJ  uJ  UJ  uj  ljuj  uj  LU  uj  uj  uJ  ui  uj  uJ 
OO  4  KOinir  4  ONON4f*«OirKN 
OOO  HO  4iAO  HOK  rim  m  H04  N  W 

r\  .h  «  k.  4>  ®  «  ho®®  «  4  «  •  4  r*  «h  ro 

noh  oin  om  ,-tin  -o  hcoiawia  o-o  or 

K®innHHc®iAroirc®iriAoo®K 

o  t  nn  «  no  ao<n  tV4  kiAiak  04  o  o 

miAK  v  k  rr>  m  n  cw«  om  «  k  4  4  c 

4  H  M  ®  H«A®  r«.H,H(V<Vr0rO4in4K 

.,  ..  .,  ..  *.  ».  ..  •  ...•■•...»•«».•  •»  •  . .  *a  •• 

oooooooooooo  OOOOOOOO 


c  o 

o  o 

LJ  UJ 

o  o 

K*  *H 

4  to 

C  00 

o  tv 
in  h 

«  V 

to  o 
o  o 


OHaCODOOOOOOOOOOCCOOOGOOO 
oo  OOOOOO  OOOC  OOCOOOOOCOOOO 

UI 


UUJUJLJUJUJUJLJUJUJ  oJU-JUJLJUJlJ^JUJUJ 
o  k  «  k  «•  H4oiA4toonnr)®o  o®  ® 
O4NNHN««KCI0  HK  OlA  4  ®  HO  4 
N  OK  NW®n  ®  n4  (►  N«®0®«  HOHO 
OOKOOK4  HNK  N  ®  O  O  IA  <0  Q  <0  K  <0  ® 

»ooo4  <o  k.  cvj  ooo  <voo«mc«KroK  m 
OOAMOO  HHHK|AOK®0®®A«A« 
®®CKC«r)<P04®<  OHHOOONtO 
OHOOOOO-CCfflCOOOOOO®®®® 

.  ..  ..  ..  ..  ..  .  ..  .,  .,  .,  .,  .,  I.  ..  .  a,  ».  .,  • 

ooo  o  o o o  o  o  o  o  o  o  o  o  o  o o  o 


in  k 
in  o 
10  c 
4  in 
0  m 

CO  -o 

O  •« 

oc 


UJ  UJ 

in  h 
tv  cv 
m  o 
O'  o 
CD  O 
CO  fK 

••i  o 
0-0 


o  o 

o  o 

UJ  UJ 

o  ro 
CM  00 
•VO* 

00  tv 
o  tv 
o  in 

•  •n 
8>  00 


OOOOOO 


>-  IT  4  K  n  f\i  HHtorc  OHr  hH  H  HN  N  N  A|  N  A'  N  N  N  N  N  A.  N  (V  f.  NNNNNNNNNNNNNAN 

h  OOOOOOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

«  I  •  4  9  l|  ■  • 

to  UJ  'UJ  u:  UJ  UI  u:  IL  LU  LJJ  UJ  UJ  UJ  LJ  UI  u;  UJ  Ui  UJ  UJ  UJ  a)  LU  UJ  LL1  UJ  UJ  UJ  UJ  UJ  UJ  U.'  LL'  LU  UJ  u:  UJ  UJ  UJ  u;  UJ  Ui  UJ  U)  LkJ  u:  UJ  LU  UI 

Z  O  H4  in  n  4K  HOt>NK®C0-flKH4NN4nnKWNOK'0HO4H4KN4HMHIAf0®IANNt> 

ui  otr  ooo  ®  riH  o<o  r-  oia,0  4  nn^  4IA  o«-oiank»k)K  4  «  a  nr  o*  oro  hoa  tv  j>  <0  ^ro  oo>  r*  n 

*-  ®  4  h  o  h  ho  oN®®iAKN44aAONcoor)OHONK  on  ho  r)  H4  ®  o  r)<on  o®  4  «-<roo*4 

.7  4  in  co  cvinr^-  hwk  roco  onninar  onnia  ia«  o*  M4  h«k)onk>ooh®iadh4h®a®  h4)  ok 

»—  K  >0  rH  IT.  K-  4  4®  Dm  4  ®  H  HOO  OOvHOtHOOO*  »H  C«J  in  K'  □  o  U  H  <0  00304344«H44490(X)00«h 

O  O  in  AO  HfOO  4  4  O  hMN  H4  4  MW  04  O  K)  ®  O®  r)  O®  ffl  IA  K  4  (OtXUO  04  14)  OPO®  O  K) 

If,  H  4  N  ro  O'  4  .HO  OK  in  HK  0«lAKOJ0inOK®4>O'0lAJ0NNNK)44444f0r)f044444WN 

K  4  CM  HA  HlA  HN  >04  K  HHNO  A  K  HHHHHHHHHHHHHHHHHHHHHHHHHH  H  H  HH 


OOOOOO  OOOOOO  OOO  OOOOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOO 


UJ  NNHHHOGCOCC  oHHHHHHHHHHHHHHr.HHHHHHHHHHHHHHHHHHHHH 

O  OOOOOOOOOOO  OOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOO 

3  •  I  11  I  II 

LLiuJUJUlUJUlUJUiUJUlUJ  UJUJUJUIUIUJUJUJUJUJUJUJLiiUJUJUJUlUJLUUJUIUJuJUJUlUJtiJUJUJuJUJUJUJUJUllllUJ 
O  re  :C  in  O  4K)41*rNcr)®N®®N4®t)ONf'JlAWn)0'®4K4lfL,»ci0'Or)fflc0ttC'0OaH44 
■y.  N  HlA  NK  4  in  *H  IT-.  HOlAJ0N0®Nin4J0ON®®40fr4O'0'0«K®>N4  0  4oKO4>(lKN0'K® 

(3  O®  N4  H  Hn  A  ®  O  N  O  A  N  N  O®  OK4  3  4)K  HN>OK  n  K  A  OJONam®lANO*HA®10*Kf90 

<  0*0®  N4  ®  HO*  N  roc  IO  CO  K  OCO  K4MM(M>K®OKK  or)N®  HK  MA«  HH4  A4  HH®  IAOK  H 

X  a  4  A  K  A  4  HAt)0*0  ®  t)A  ®0»A  OCK  4  r\IO>  OlA  HO  JO  IA  □AO>NHHo®r)Nn®OHHNO"0 

AAflttfBWAN  »f  O  O  HHAX.aOON  H  <0  HK  040  r!  C  W  'C  'C  t  r?  N  ff  4  J  CO  4  C  NIOC  C  O 
K4AlA>0r0l0  4AlA«  O  OfO  OO  K)  K  N  O  O  HHHOHOOKAAA«K®®KKO<IOKK®®K  <Oin 
N®  HWK  HNW41A®  C  HHHHNNWWW444444KJWWOWr  «  W  K)  W  W  W  H  W  IO  «W  M  W  M  O 


OOOO  O  O  OO  OOOO  OOOOOOOOOOOO  OOOOOOOOOOOOOOOOOOOOOOO  oo 


HHHHHHHHHHHH HHHHHHHHHHHHHHHHHHHHHHHHOOaO  OOqOhHHH 
OOOOOO  OOOOOO  OOOOOOOOOOOO  OOOOO  OOOOOOOOOOO  OOOOOOOO 


X 


UIUJUJUJUJUJUJUIUIUJUJUJ  uJUJUJlJUJUJUJUJUJUJuJliJUJlillUlilUKJIUJUJiJIUlUlIttUJUJUJUJ 
OOOC  OO  OOOOOO  OOOCOOOOOOOOOOOOOOOOOOOOOOOO 
OOOC  OO  OOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOCOOO 
coco  co  ro  c  c  oc  occoocrocoococeocccccooccocc 
n a oo  o o on oo  no ooonoooon n no ooooo onnnoonooQo 
cm  C4  c  ir  oin  om  cur  oiroir  om  or  oo  cm  omom  cm  coco  cmcooo 

ok  an  ok  in  rv  ok  in  tv  ok  in  tv  on  an  ok  an  okanokanok  anqa  om 

CN  Ar  H®  «4  NO  KIA  K,  C®  0  4  HOP  A  A  C«.  «»T  HO  K  4  N  CC  If  K5  HOP  intV 

OOCOffl®  \NNN4)<0  4)  4)  A  A  A  A  A  4  4  4  44K)K)Mr)NNNNNHHHHO‘®4« 


UllUUlUJ  lilUJ  UJUJ 
OOOOOOOO 

o  o  O  o  o  o  o  o 
O  CoC  oo  oc 
noooooon 
c  c  o  c  it  c  r  c 

in  om  o  tvm  n  o 

non  c  .h  r:  in  «r  c 

N  4  OO 


C  C  OOOO  OO  O  O  OO  ooo  OOOOOO  O  0900090  9  0  9  0  900090C-»0  9  0  0  9  0  C  99 

n  it  it  it  it  ii  i  !»•(••  it*  il  •  il- »  •  i|  ill  ii -11*  it*  it  <  »•  n  it  >il  >il«  ii  il  :l  umi  ii'it'U’il'ii'Uit  «J  it  i  it  mi  -ii 


F* 


HN*OT  AC  N®  »  O  HN  K)4A  CKOC  OHNr)4  m  -OKCQC  OHNr)4lA®K®0>  o  HN  tO  4  IA  ®K  B  »  O 

H  HH  H»4«iH«lH«tNNNNNNNNNNnnnnnnnnnn444  4  4  4  4  4  4  4  m 


148 

I 


25J000E  01  0.3535AS64E  31  o;i25o2H8oE  0?  0-84797474E  00  0-7i9o6ll6E  00  0.34667 


I 


I 

I 


□ooaaooooooo oooooooa ooaQQoOODrt h 

a  ooaoaaoooao  oooooooooooooooooo  a 

UiuJUULJaiUUULJiJUULJUULJUJUiuiULJLiLJULiUiuUJiJ  LJ 

^wwwrxNjrfoHcoiDotHK)^  on  ory  rv  fv  k>  *©  mo  o  o 

<  -O  ON  «3  MlA  o  HO  H  Ht  OOHHH8)  HK)  NlH  V  «  HOOo  o 

-J  .'X  O'  "0  C  HO  L  HN  NNn  »  N  no  O  9>  fi  90  o  a 

a  (vn  H(v«cvj-otf,  rt  vHtHotf.O"C»c  o'o^ojro^oocc 

<  oos  innK)  COW  HOOMI(V«ONVN»9(M>(M>00 

N*  fs  W  H  HO  orJNC  ON  COiOlT.  at  <ON  ISO  C  O'OOOOC'GO 

•  rv  c  H  «  r«0  H  <C  HJN  O  (Vi|n4N<00>0lOO  OOOOOOOOOCD  o 

X  f«UMO<H;NMCD(000(MM>0>0  >  O  O  O  O  O  >  O  O  O  0>  H  H 

O'  ■  -  •  •  * . • 

C  o  c  a  o  o  c  o  o □ o  o  o  ooocoaoc  coooooooac  o 


I 

I 


I 

a 


a 

5 


H  H  H  H  H  (V  HHO  O  oo  ooo  HOOOHHHHH  H  H  H  H  HH  O 
oaccocooooac ooc oooooocoooooooo c 
I  (I  U  I  i»  ll  I*  >4  il  il 

uJ^iiLiJUJLL:aJuiUJU.a.UiLLJLUuii*JUJUJUJUiu;uJUJUJu.'Ulu.'UiUiUJuiuJ 
n-  cn^cc  rv  m  rv  rx  cv  cc  r-<  m  «  c  □  n  n  t  i\j  t  itcd  ^  in  c*  »-» ®  «-*  *4®  «o 

l/>  U«  y  on^fVJt  «CON  *  O  V  X*  OQON^O  ON  (O  t  O  O  OK)  CO 

«  urjaj  rv  a.  mn  v  *h  r-  rv  oau  mb  Ni  or)  a  m  o  w*o  ho*  hid  in 

o  o  ao  offi  <r  oottUAriffao  cooo©  rotv.  on  ©o  »  o  air.  fx 

(m  a  oj  aa  v  a  omNo  v  o  ra  ©  ry  cc  ©  w  aa  ©  ry  n  ax  t-  o  njo>  rx 

at>  h  o«o  «  ot\i  own  ©noo  ot>  oo  on  in®  mn  to  imns  fx 

wyro  ®  C  C  HC  «  (VO  «»,lMNM£r»INKOCac»IKlN©  rv. 

rx  o  (N  f.  N)  K)  (4<o  »K4  rj  n  kj  N)  fa  u  nro  r-  »i  h  r j  n  ra  h  ©  cj  cith 


coioaoooocooc  cacocoaooooaoaococ  a 

iiiiiiiiiii  ill  I  I  I 


DcaOCCCCOOOOGOC  OrlHrlrlHHtMMMWt  *  m  O 
o  o  o  o  o  a  ooooooooaoaoooooooooooao 


r 

cr 

o 


N  K)  O  K>  fO  rx  inri  HH  v  HN  H  O  OIA  €  H©  ID  or  HO©  DH 

If .  (\j  ©  (D  flO  Hf  O  If  OOtMMMMO  H«  K)  HONM« 

in  tNin  iMwvo  Hto  oroxn  H«in«t  oh»  v  ©>▼«  o« 
'O  hho  roo  in  o  nnuin  m®  oot  hoia  nn<M«  o  hnhmw 
mo  rv  *  c  v  oao  orjif.  ©  ho>  ry  ©  ©  «h®  mc  o  m  rv  kj  ®  a  ®  cv 

H  fON  CNlA  HO  A4DON  to  OOlft  (MOfiN  O  04  a  O  O'  HO  If, 

▼  o  ©>  ■»  ©insoNoto  t  Hcvm  »inr)N  n  h®  HHHHt  wn 

N*  NMOOOOOO®NX-tM(\JHO©tWHH«K)  HfON  rt  Nt 


oaoooooooooo  ooaooooaoooooooooo  o 


oooooaocaooaococooaoooHHHH(V(\i  t\ir< 
coocoaaoaooccjoaocaooooaoocoooo 

LJuiOiUiUlOiLJ^lUJUJUlUJaJLjm^UJUJUJUJUJU/UiLJlDUJUJUJUiUi 
c  t  w©  oir  on  ©  n  o  on  t  hk)  o(V  h  om  h©  oroo©  CO 

•  «©(VHfQnjinONHDW(\JtOH©^HM>HON»OONON 
13  lTB©©WaNONIOt«flKf*et«r©OH«B©NffHr 
*<  (vhtnoono'Dcon  o©omnin«  H©inHinomoHON 
»  o  o  rv  rx  ©t  o»  nw  on  «o  oiMnacinniMroHHWircN 

hhhimon  mo©  Hino  o  v  on  mmN  ^  oj®®o  cv  ®  ®  in  -vo 

•  •*  i>  r)  ®  m  »x  «  »  «  w  «  «®  o  «  in  to  nr  o  ®  r>  o  «m  r>  #®  r*  *•  ►  w*® 
x  x  Doooo>oo>»o«N©invniON(MHHHomr)H©nH« 

CE  •.  •.  •.  •.  •.  i.  •.  •  «,  •,  •.  t.  •.  •.  •  •• 

O  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


>-  MV  (\  VIOl  A  (V  (V  O  (V  (M  (V  HHrtHi-rOOnOCiHHtyNrtin 

►*  c  oooooooooooocoooaoooaooooacoo 

*-  II I  III  II  •  II  •  II  ll 

tr.  ujuiujuj  LUttiaJLUir'UijUUJiiiuiiiiiiiU’iiJaju:  uiujlull:  uuujujujujuj 

7  ®rx  H(\J  t  N  MO  t  HN  «  Ht  CD  ©  (  m  N  n  N  r\j  IT*  ®  O  N  W©  IMIO 

u;  o  tN  n  or-  rytccoirtofjNt  o  t  cn  ©  on  ®  <n  r*  k: 

►-  DNry®Honoo>o(\j©inDttNMnin©(y©o 

zz  ok)®  Ht  no  ©minmf'i©N  orn  *h  *h  o®  kin  n  hn  in  (m»d  imt 

O  ONin  (NlHO  ffl  HOHO  OO  O  *H  *H®  V  IT'  <0  ®  V  ^  rxi/>  »H  »H  ®  N. 

««OK)0©OK)0*t  OKJiotf  HK)K)HOf  tOK)  HCMnm  ©© 
NKMn©©©©N©mnoNm©©NHmN  ooh«  ownit  ho 
HHHHHHHHHHHHNiAfONHHNV  ION  Hlf>  *H  lT»  HN  tN 


I 


0  00000000000000000000000000000(1 


111  HHHHHHHHHHHrl  HHHHHHOOOC  OO  O  H  H  H  N  C\1 
o  ooooannoonooooonaooooooDODODoo 
n  .i  ii  ii  t»  i 

►*  UJUI*llU«L!ja4UILliUldiUlUJliJUIllJaJliiUJLliUIUIUJUJUiUiUJUlUlUlUi 
—  ©  tc«cr  K)n  ©  HH®  ON  t  rl©‘0  N©N  HON  rtfx  OC  H  NlT  NK) 

z  ©^tt>otffl©Non»mN©©NK)toHmNmvfflNin^N 

O  N  nN  >ONHN©OntN  O©  O  NN©  ©  N(ND  HHf  NNO 

<  NIOONNC'  OONO>(MVNVMOON<Or)<MONO«HCD'r(M(Oa3» 

x  a  m  k)  «o  rtin  o>(\j?  n>  on  (N«in©8>oo©iAH^inNin«o 

c>  hc  t  on  h©  riryo  com  m  r  i  ri  fl  ©  m  Hfy  in  n  ®  t  n  o  inm 

lA  N®  o  HO  HH  H»  ©  NNK)©>  ©  K)  0©C  Ifl  iA  t  K)  KinifMA  «N 

•O  K)K)  ^  *  t  ttKirOKJNNHH  HH®  ®  U>  ^  K)  CV  -*rx  K)  ®  CM 


OOOO  OOOOOOOOOOO  ooooooooo  oooooo  o 


a 

a: 


o  O  O  O  OOOOOOOOOOOOOOOOOOOOOOOOOC  z> 


I 

1 

I 


UIUI4UUI  UIUI  111  Ul  m  HI  HI  LAI  UI'H  U1  HI  ym  111  Hi  Hi  HI  114  111  UlUlllllUUiUl  ill 

OOOOOOOOOOOOOOOOOOOOOOOOOOOO.  oo  o 
OOOO  OOOOOOO  OOOOOOOOOOOOOOOOCN  oo  o 
O  C  CO  OC  OCOOOOCCOCOOOC  O  O  O  C1  o  c  ot>  ooo 
oooo ooooopoonoooooooonoooooo ooo 

cm  c®  oir  oir  oncir  cr  ononDir  or  or  oin  ct  oir  c 

>AN  ONiON  ON®  N  ON®N  ON®N  ON®  N  ON  ®  N  ON®N  O 
N  t  f-  O  HK-©CCMnNOH»<CCr-lfNON«©f  n«®N  C. 

NNNNnnnnt^V^t®®®®©©©©«NN  .x  N  «  ®  ®  ®  o 

000000000000900003000000000300 9 


Ul 

a 


HN«  ♦  ®  ©  N  ®  9  OHNKIt  IT  ©NO©>  OHN^  ©  ®«N«90  H 

®®®®'n®®®®©©©©©©©©«©NNNNNNNNNN«© 


z 


149 


r> 


1 


u 

2 


a 
c ) 


G. 

a 


2  H 

Z3  er 

o 


4  t4UJ 

♦loo 
I  o 
♦  I  UJ  o 

•  h  M  o 
I  «  o 

•  *oo> 

•<©  . 

1  ro  u 

1  o 

•  •  *0 

»  r.  -4 

I  •  c 

•  o 

I  M 
»  ^ 

I  CM 

•h  rs. 

•  (V 

•  N 

1  ra 

•  ON 

•  in 

•  ui  • 

•  H  O 

»  O 

I  O 
HO1 
I  O 

•  o> 

1  O 

•  ♦ 

I 
I 
I 


© 

IT. 


© 


ff 

H 

• 

V*  <pH  1 

■f*H  H  r  II  Vt 

1 

1 

1 

1 

a  in 

ON 

U.  Z 

O 

01 

01 

[0 

01 

to  0 

ui  ui 
O  CM 

UJ 

M 

UI 

UJ 

►  a 

-i 

3b. 

r>  • 

N-  © 
•-•N 

O* 

© 

O 

*• 

H 

CO 

© 

© 

if.  O 

UJ 

a  a 

UJ 

.  •* 

2  O' 

0  <o 

•< 

X  © 

•  v 

CM 

© 

rs. 

w 

© 

n 

n 

n 

« 

© 

H 

O 

CM 

© 

H 

O 

© 

CM 

u.  x: 

• 

O 

0 

0 

O 

O 

• 

I 

I  OLA 
«  Q« 
I  CM 

•  UI  • 

I  CM  O 

1  to 

•  H 

«H© 

•  H 
I  05 
I  H 

1  ©  O 

I  •  O 

•  o 

I  ill 

1  fc 

•  00 

*"«  *H 

I  © 

•  rl 
I  HC 

•  O  H 

«  © 

I  UJ  • 

•  ©  o 
I  ©  I 

»  w 
•H© 

I  V 

•  © 

t  o 

•  Cl  «H 
I  •  c 
1  o 

•  OLil 

•  H 

•  o 


•  HO 
»  O  C 

•  » 
I  ttJ-  • 

1  II  N  O 

•  CM  • 

•  N 

^  CM 

I  N 
I  CM 
I  N 
I  i/>  H 
I  •  O 
I  =» 

•  IU.' 

1  o 

I  © 
*H  W 

•  © 

•  n 

•  H« 

•  OK 

♦  •  N 

♦  I  UJ  • 

♦  lac 

•  on 

♦  •  o 


Z2 


1  HH 

«4  «H  | 

1 

1 

1 

m 

•IHHHil-HH 

♦  HO 

•4  Hl  »  HH  1  O 

O 

01 

01 

00 

s  s? 

0  s 

UI 

H 

»o 

© 

© 

« 

O 

(M 

UI 

© 

© 

O 

© 

© 

« 

H 
•  m 

UI 

O 

0» 

N 

Oh 

O 

© 

CM 

vf 

ill 

© 

CM 

w 

© 

0 

n 

3 

0 

UI  f, 

CM 

• 

« 

fs. 

O 

© 

H 

* 

O 

O 

O 

•  N 

O 

•  *  .  •  * 

O  O 

150 


•  i*  rs  o 

•  I  (M  • 

I  N 
w4  tM 

i  r^ 

•  ru 

i 

i  r<H 
l  •  o 
I  a 
I  l  Ul 
»  « 

♦  I  c 

♦  *4  to 

•  « 

♦  i  to 

♦  I  H-O 

4  an 

♦  i  ui  • 

♦  •  b  oa 
i  a  ti 

♦  l  a 


v*  •  *4  v4  V  O 


n'.f.  results  for  round  trip  „  . 

NUMBER  OP  POINT*:*  »1  X  INCREMENT*  0.16363636E  0*  T  INCREMENT*  n .  20C  (1 0  o  0  CE- 0 1 


r 


♦ 

♦  ‘ 


O 

I  o> 


rtU. 

O  o 


♦  O 

<0  O 

no* 
«  • 
n  o 
o 
n 

N  H 


cm 

fv 

cvi 

r-«* 

r<CM 

ON 

4fl 

(1/  • 


If* 

IT\ 


om 

o» 

CM 


m 

c  o 


•0® 
O  H 

«o 


ir\  o 
m  • 


nr 

‘  * 

in 


o 

(I  Ul 


o 

O* 


H® 

o  o 


OJ-  • 
rn  o 
CM  II 
N 
wit* 

h 

CM 


fO 

o 

f-l 


Ul*  • 

O  O 
O  II 
o 


a  HH  HH  »  ******  **  •  •  *******  •  «•  •  ******  9  *4*  *** 

^CCOCDOC  Q. 

O  O  O  O  oooo 


0  *l» 

•  MH-v*  *«  **'  9 


13  t|| 
<  O 
X  P» 


X  o 
a c  9 
o  o 


Ul 

o 


3 

£ 


152 


I  HUJ 

I  o  o 
I  o 


♦ 


♦  •  f»  o 

♦  I  O  CD 

I  K»0» 
«  •««  - 

♦  INO 
I  « 

>  I  w 

I  fN  H 
I  •  «3 


r 


0 

ro 

•c 


z 

ui 

*: 


a 

a 


a 

Z  H 
3  «C 
o 


O  O 
O 
UI  o 
V  o 

<0  o 
wo» 

H<0  • 

w  a 

0 

w 

r-  h 


CM 

•* 

(M 

N 

HM 

or* 

0 

Ui  • 
H  t) 

o 


sD 

a 

r* 

o 


O 

o> 


a: 

u 

2 


0 

0 

ft 

1 1\ 

ft 

om 
a« 
CM 
UI  • 
CM  O 
0 


to 

H 
« 
•4 
H«0 
o  H 
0 
Ui-  • 
0  O 
0  t 
ft 
0 
* 

0 

ft 

CM  *-* 
•  O 

o 

II  UI 
w* 

o 
o 
ft 
o 
H» 
O  CD 


CM  • 

K 

CM 

fs 

CM 

r- 

0  T* 
•  0 


to 

0 

ro 


I  • 
»o 

I  II 


•4  v«*  •  V*  ;«4  I  -w4  <v4  wi  v*  I 


a 

N 

• 

1 

i 

t 

• 

«r4vt  vt««  1  v« 

l 

* 

i 

t 

fl  1  *l»t 

•4  •  <•*  »4  ;w4  wi 

t  • 

a  tr 

H 

rl 

*4 

CD 

O 

CD 

C. 

H 

(4 

o  r 

UI  9 

O 

CD 

o 

o 

O 

o 

O 

o 

o 

u.  z 

0 

m-m 

ft  UI 

ai 

m 

UI 

Ui 

Ui 

UI 

lit 

ul 

ui 

to  o 

x  ft 

o* 

0 

CM 

0 

0 

0 

o 

•0 

►-  u 

a  o 

0 

90 

0 

»h. 

CM 

CM 

H 

ft 

P»k 

— ; 

90 

o 

r*. 

o 

r* 

0 

0 

ft 

CM 

uu. 

-Of 

r* 

ft 

o 

CM 

0 

ft 

0 

yv 

to  a 

X  » 

0 

CM 

0 

H 

ft 

O' 

0 

ft 

CM 

ui 

ft  r* 

0 

ft 

0 

0 

0 

H 

CM 

0 

H 

a  cr 

o  n 

0 

CM 

ft 

O 

90 

0 

ft 

ft 

Hx 

ui 

Z  CM 

wt 

H 

CM 

90 

0 

H 

0 

CM 

-  m 

»  • 

.  » 

u*  X 

o 

o 

O 

O 

O 

O 

O 

O 

0 

g 

•  *3 

•  9 

• 

•  9 

•  • 

• 

zz 

i  i 


[  I 


154 


o 


2  FOR  CASE  8,45 


starting  pass  iuMMnfcB 

AR£A(S)  AT  START  OF  MSS 

2D  ARM  3.1*»16770E  03 

AFTER  I  NTFRP0lAt  I  ON 

2D  ARM  3.1W4971F  03 

ATTER  OUTPUT  mJBRpR 

2D  *8e*  0.7R7457RE  02 

AFTfR  Sf6p.BY*STtP  PROPAGATION 
2°  *Hp4  0^<»7lft69t  02 

AT  ffc  p  mE0!  up 

2D  ARM  0,320«774F  03 

ATTER  STtP»BY.STfcp  PROPAGATION 
2D  ARM  3  : 3209204F  03 

AF  TER  BACK  Hi  PROP 

2D  A»f-A  0  ;  2973o94fc  oi 

AFTER  SYEP.BTsSTLp  PROPAGATION 
2D  area  3i2V7309Bt  03 

AFTER  MEDIUM 

2D  AREA  0.670*>350E  03 

after  step-by-step  PPOTAGATI  0'.' 

2D  AREA  C.6705364E  03 

finished  pass  number  2 

POWER 

NEAR  FIELD  POWER*  0.S705364E  03 

OUTPUT  P0WEO  after  MIRROR.  0.3968960E  03 

PDwE«  oIFFebENCB*  0.47363U4E  03 
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APPENDIX  E 


LISTINGS  AND  FLOWCHARTS  SSCP(ECS) 


0-2f85- H 


^  PR51  ^ 

^ - 1 - J 


Fig.  El.  P R 5 1  logic  flow  diagram. 
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waevmm 


(expansion  of  PR51 
round  trip  pass 
block  of  Figure  El ) 


P50INT 
reduce  wave(s) 
to  compact  form 
to  begin  pass 


reflect 
wave(s)  from 
back  mirror 


reflect 
wave(s)  from 
output 


X.  back 

medium  mi  rror 

location?  ^  1 


elsewhere 


id  lie*' 
lutjttnn 


PSOSBS 
P root  gate 
to 

medium 


e  1  sw where 


edi  vi 

exists? 


(return  to  Figure  El, 
PR51  logic 
flow  diagram) 


Fig.  E2.  Logic  flow  diagram  for  round  trip  section  of  PR51 
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APPENDIX  F 


GLOSSARY  OF  PRIMARY  SSCP(ECS)  NAMES 


ALLOUT 
AMPINT 
AREA21 
A  SR 

CASE 

CASENO 

CAVITY 

CSCALE 


El  A 1 
DIA2 
DXOLD 

GAIN 

GAINA 

GZERO 

IFILE 

INFOUT 

IRC  ODE 


Subroutine  by  which  GPOP  output  is  produced 

Array  for  defining  initial  amplitude  of  waves 

Subroutine  for  computing  area  of  a  wave 

Input  variable  defining  pseudo-three-dimensional 
parameter 

Input  variable  used  for  titling  output 

Namelist  block  name  for  describing  number  of  cases  to 
be  analyzed 

Namelist  block  name  for  printing  cavity  description 

Subroutine  for  determining  printer  plot  scale 
increments 

Input  variable  defining  diameter  of  output  mirror 

Input  variable  defining  diameter  of  back  mirror 

Array  saving  all  DX  values  for  each  step  of  the  last 
half-pass  through  the  medium 

Array  of  gain  and  refractive  index  variations 

Subroutine  defining  GAIN  array  for  medium  propagation 

Input  variable  initalizing  gQ  for  GAINA 

Input  variable  initializing  normal  output  file  code 

Array  defining  output  option  for  near  field  printout 

Input  variable  describing  file  code  from  which  restart 
information  is  to  be  read 
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IROUT 

IRTAG 

LBLPOS 

MEDIUM 

MIRROR 

MIRRRS 

NCASES 

NORM2D 

NPASS 

NWAVES 

NWDIM 

NX 

NXC1 

NXC2 

NXCOUT 

NXDIM 

NZDIM 

ONEOUT 

OUT1-OUT7 


Input  variable  describing  file  code  to  which  restart 
information  is  to  be  written 

Input  variable  describing  pass  number  of  wave(s)  to  be 
taken  as  restart  information 

Subroutine  positioning  column  heading  label  for  GPOP 

Namelist  block  name  for  printing  medium  description 

Subroutine  for  reflecting  waves  from  mirror 

Namelist  block  name  for  printing  mirror  descriptions 

Input  variable  setting  number  of  cases  t  be  analyzed 

Subroutine  for  preserving  2D  Area 

Input  variable  setting  number  of  passes  for  a  case 

Input  parameters  defining  the  number  of  propagating 
waves 

Dimensioning  parameter  for  maximum  number  of  waves 
that  can  be  handled  by  the  SSCP(ECS)  or  SSCP 

Input  parameter  defining  the  number  of  mesh  points 
along  the  x-direction 

Computed  parameter  defining  the  number  of  x-intervals 
across  the  output  mirror  at  start  of  pass 

Computed  parameter  defining  the  number  of  x-intervals 
across  the  back  mirror 

Computed  parameter  defining  the  number  of  x-intervals 
across  the  output  mirror  at  end  of  pass 

Dimensioning  parameter  for  maximum  number  of  mesh 
points  along  x  that  can  be  handled  by  the  FICP 

Dimensioning  parameter  for  maximum  number  of  last- 
pass  intensity  stores  possible 

Subroutine  for  defining  vector  from  GPOP  options 
Arrays  used  by  GPOP  for  defining  values  to  be  produced 
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P50GAI 

P50GET 

P50INT 

P50MED 

P5CSBS 

PGPLT2 

PR51 

PWROUT 

RAT3DA 

RCURV1 

RCURV2 

REF  LI 

REFL2 

RSTART 

RUNOPT 

SETMIR 

SO 

SONEW 

SOTOT 


Subroutine  for  computing  set  of  gain  values  at  a 
z-location 

Subroutine  for  retrieving  a  value  from  the  last  pass 
intensities  matrix  (SO) 

Subroutine  for  reducing  the  expanded  waveform(s)  to 
begin  ?.nother  pass 

Subroutine  for  propagating  one  way  through  a  medium 

Subroutine  for  propagating  through  an  empty  cavity  (no 
medium)  region 

Subroutine  producing  GPOP  printer  plot 
Main  program  of  SSCP(ECS) 

Subroutine  computing  near  field  and  output  power 
Subroutine  for  preserving  3D  Area 

Input  parameter  defining  the  radius  of  curvature  for  the 
output  mirror 

Input  parameter  defining  the  radius  of  curvature  of  the 
back  mirror 

Input  parameter  defining  the  power  reflectivity  of  the 
output  mirror 

Input  parameter  defining  the  power  reflectivity  of  the 
back  mirror 

Subroutine  forming  an  input  wave  by  restarting  a  previ¬ 
ous  run 

Namelist  block  name  for  printing  run  options 

Subroutine  for  setting  up  a  mirror  vector 

Array  storing  last-pass  intensity  values 

Array  storing  current  intensities  until  they  can  be  stored 

Array  storing  sum  of  old  and  new  intensities 
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T 

TILT  1 

TILT2 

TLOSS1 

U 

V 

VR1 

VR2 

VALUE 

WAVES 

W  OMEGA 
X 

XCENTR 

XIZERO 

XKZERO 

XLAMDA 

XMIR1 

XMIR2 

XWIDTH 


Array  (dummy)  used  in  propagation  routines 

Input  parameter  defining  mirror  tilt  angle  of  output 
mirror 

Input  parameter  defining  mirror  tilt  angle  of  back 
mirror 

Input  parameter  defining  power  loss  of  output  mirror 

Array  storing  propagating  waveforms 

Array  (dummy)  used  in  propagation  routines 

Virtual  image  source  radius  of  curvature  for  output 
mirror 

Virtual  image  source  radius  of  curvature  for  back 
mirror 

Subroutine  computing  one  of  several  values  for  complex 
numbers 

Namelist  block  name  for  printing  propagating  wave 
characteristics 

Array  defining  gaussian  w  for  each  input  wave 

Array  defining  physical  position  of  each  field  point 

Input  parameter  defining  the  x-value  for  the  center  of 
the  propagation  array 

Input  parameter  defining  I  for  GAINA 

Array  storing  Kq,  or  2tt/\,  for  each  wave 

Array  defining  X.  for  each  wave 

Array  representing  output  mirror 

Array  representing  back  mirror 

Input  parameter  defining  minimum  spatial  extent  (in  x) 
of  the  beam 
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ZLEN 

ZLENMD 

ZMDLOC 


Input  parameter  defining  distance  between  the  output  and 
back  mirrors 

Input  parameter  defining  the  length  (along  z)  of  the  gain 
medium 

Input  variable  describing  position  of  medium  in  cavity 
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